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Abstract. We develop a model of an electrorheological fluid such that the fluid is con- 
sidered as an anisotropic one with the viscosity depending on the second invariant of the 
rate of strain tensor, on the module of the vector of electric field strength, and on the angle 
between the vectors of velocity and electric field. We study general problems on the flow 
of such fluids at nonhomogeneous mixed boundary conditions, wherein values of velocities 
and surface forces are given on different parts of the boundary. We consider the cases 
where the viscosity function is continuous and singular, equal to infinity, when the second 
invariant of the rate of strain tensor is equal to zero. In the second case the problem is 
reduced to a variational inequality. By using the methods of a fixed point, monotonicity, 
and compactness, we prove existence results for the problems under consideration. Some 
efficient methods for numerical solution of the problems are examined. 



1. Introduction 

Electrorheological fluids are smart materials which are concentrated suspensions of polar- 
izable particles in a nonconducting dielectric liquid. In moderately large electric fields, the 
particles form chains along the field lines, and these chains then aggregate to form columns 
(see Fig. 1 , taken from [18]). These chainlike and columnar structures cause dramatic 
changes in the rheological properties of the suspensions. The fluids become anisotropic, the 
apparent viscosity (the resistance to flow) in the direction orthogonal to the direction of 
electric field abruptly increases, while the apparent viscosity in the direction of the electric 
field changes not so drastically. 

The chainlike and columnar structures are destroyed under the action of large stresses, 
and then the apparent viscosity of the fluid decreases and the fluid becomes less anisotropic. 

Constitutive relations for electrorheological fluids in which the stress tensor a is an 
isotropic function of the vector of electric field strength E and the rate of strain tensor 
e were derived in [19], and for an incompressible fluid there was obtained the following 
equation: 

a = -ph + a 2 E <g> E + a 3 e + a 4 e 2 + a 5 (sE <g> E + E <g> eE) + a 6 (e 2 E <g> E + E ® e 2 E). 

(1.1) 

Here p is the pressure, I\ the unit tensor, on are scalar functions of six invariants of the 
tensors e, E ® E, and mixed tensors; are to be determined by experiments. 

In the condition of simple shear flow, when the vectors of velocity v and electric field E 
are orthogonal and E is in the plane of flow, the terms with coefficients at2, «4, as, «6 g iVe 
rise to two normal stresses differences (see [19]). But these terms lead to incorrectness of the 
boundary value problems for the constitutive equation (1.1), and very restrictive conditions 
should be imposed on the coefficients «2, «4, «5, a$ in order to get an operator satisfying the 
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conditions of coerciveness and monotonicity (the condition of coerciveness is almost similar 
to the Clausius-Duhem inequality following from the second law of thermodynamics, and 
the condition of monotonicity denotes that stresses increase as the rate of strains increase) . 

The constitutive equation (1.1) does not describe anisotropy of the fluid; in the case of 
simple shear flow (1.1) gives the same values of the shear stresses in the cases, when the 
vectors of velocity and electric field are orthogonal and parallel (a is an isotropic function 
of E and e in (1.1)). 

Stationary and nonstationary mathematical problems for the special case of (1.1) are 
studied in [20]. It is supposed in [20] that velocities are equal to zero everywhere on the 
boundary and the stress tensor is given by 

a = -ph + 7 i((l + |e| 2 )^ - 1)E ® E 

+ (72 + 73|£| 2 )(1 + \e\ 2 )^ e + 74 (1 + |e| 2 )^ (eE ® E + E ® eE), (1.2) 

where |e| 2 = Y17j=i £ iji n being the dimension of a domain of flow, 71 — 74, are constants, 
and k is a function of \E\ 2 . 

The constants 71 — 74 and the function k are determined by the approximation of flow 
curves which are obtained experimentally for different values of the vector of electric field E 
(see Subsection 2.2). But the conditions of coerciveness and monotonicity of the operator 
— div(cr + pl\) impose severe constraints on the constants 71 — 74 and on the function k, 
see [20] , such that with these restrictions one cannot obtain a good approximation of a flow 
curve, to say nothing of approximation of a set of flow curves corresponding to different 
values of E. 




Figure 1 . Fibrous structure formed by the electric field for alumina particles 

Below in Section 2, we develop a constitutive equation of electrorheological fluids such that 
a fluid is considered as a viscous one with the viscosity depending on the second invariant 
of the rate of strain tensor, on the module of the vector of electric field strength, and on the 
angle between the vectors of velocity and electric field strength. This constitutive equation 
describes the main peculiarities of electrorheological fluids, and it can be identified so that a 
set of flow curves corresponding to different values of E is approximated with a high degree 
of accuracy, and it leads to correct mathematical problems. 

In Section 3, we present auxiliary results, and in Sections 4-8 we study problems on 
stationary flow of such fluids at nonhomogeneous mixed boundary conditions. Here we 
prescribe values of velocities and surface forces on different parts of the boundary and 
ignore the inertial forces. The cases where the viscosity function is continuous and singular, 
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equal to infinity, when the second invariant of the rate of strain tensor is equal to zero, are 
studied. In the second case the problem is reduced to a variational inequality. 

By using the methods of a fixed point, monotonicity, and compactness we prove existence 
results for the regular and singular viscosity functions. In the second case existence results 
are obtained at more restrictive assumptions. Here the singular viscosity is approximated 
by a continuous bounded one with a parameter of regularization, and a solution of the 
variational inequality is obtained as a limit of the solutions of regularized problems. 

Section 9 is concerned with numerical solution of the problems on stationary flows of 
electrorheological fluids with regular viscosity function. We consider here methods of the 
augmented Lagrangian, Birger-Kachanov, contraction and gradient. 

In Sections 10 and 11 we study problems on flow of electrorheological fluids in which 
inertial forces are taken into account. Here we consider nonhomogeneous boundary condi- 
tions in the case that velocities are given on the whole of the boundary and in the case that 
velocities and surface forces are prescribed on different parts of the boundary. With some 
suppositions existence results are proved. 

2. Constitutive equation. 

2.1. The form of the constitutive equation. It has been found experimentally that the 
shear stress and accordingly the viscosity of electrorheological fluids depend on the shear 
rate, the module of the vector of electric field strength, and the angle between the vectors of 
fluid velocity and electric fields strength (see [18, 22]). Thus, on the basis of experimental 
results we introduce the following constitutive equation 

aij(p,u,E) = -pSij + 2<p(I(u), \E\,fi(u,E))£ij(u), i, j = 1,... ,n, n = 2 or 3. 

(2.1) 

Here, <Jij(p,u, E) are the components of the stress tensor which depend on the pressure p, 
the velocity vector u = (u±, . . . , u n ) and the electric field strength E = (E±, . . . , E n ), Sy is 
the Kronecker delta, and eij (u) are the components of the rate of strain tensor 

Moreover, I(u) is the second invariant of the rate of strain tensor 

n 

I(u)= 5>y(u)) 2 , (2.3) 
and if the viscosity function depending on I(u), \E\ and fi(u,E), where 

So fi(u, E) is the square of the scalar product of the unit vectors ^ and |Jjp The function fi 
is defined by (2.4) in the case of an immovable frame of reference. If the frame of reference 
moves uniformly with a constant velocity u = (u\, . . . ,u n ), then we set: 

d^rir (»> 

As the scalar product of two vectors is independent of the frame of reference, the constitutive 
equation (2.1) is invariant with respect to the group of Galilei transformations of the frame 
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of reference that are represented as a product of time- independent translations, rotations 
and uniform motions. 

It is obvious that fi(u, E)(x) £ [0, 1], and for fixed y±, yi £ R+, where R+ = {z £ R, z > 
0}, the function y% — ► 1/2,2/3) reaches its maximum at 2/3 = and its minimum at 

?/3 = l when the vectors u(x) + -u and E 1 are correspondingly orthogonal and parallel. 

The function fi defined by (2.4), (2.5) is not specified at E = and at u = 0, and there 
does not exist an extension by continuity to the values of u = and E = 0. However, at 
E = there is no influence of the electric field. Therefore, 

<p(yi,0,y3) = <p{yi), ys € [0,1], (2.6) 

and the function fj,(u, E) need not be specified at E = 0. Likewise, in case that the measure 
of the set of points x at which u(x) = is zero, the function /i need not also be specified at 
u = 0. But in the general case we should specify fi for all values of u. Because of this we 
assume that the function /x is defined as follows: 

„ w . / al + u(x) + u E(x) \ 2 

u(u,E)(x) = { — = — r-V^ r, (, ) , (2.7) 

PV A 7 \ay/n + \u(x)+u\' \E(x)\ v ; 

where / denotes a vector with components equal to one, and a is a small positive constant. 
If u(x) / almost everywhere in Q, we may choose a = 0. 

2.2. Assumptions on the viscosity function. Flow curves of electrorheological fluids 
obtained experimentally for /j,(u,E) = have the form as displayed in Fig. 2 (cf.,e.g.,[22]). 
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Figure 2 



These curves define the relationship between the shear stress r = a±2 and the shear rate 
7 = £i2(u) = \ ^ for a flow that is close to simple shear flow. Line 1 is the flow curve for 
\E\ = 0, and lines 2-4 represent the flow curve for increasing \E\. 

Flow curves are obtained in some region, say 70 < 7 < 71, 70 > 0. Experimental results 
for small 7 are not precise, and one has to extend the flow curves to R + . It is customary 
to extend flow curves by straight lines over the region 71 < 7 < 00. One can prolong flow 
curves in [0,70) such that either r = to for 7 = or r = for 7 = (see the dash and 
dot-dash lines in Fig.2). 



The viscosity 7/(7, E) of the fluid is determined as 

ri(l,E) = ^, (2.8) 

and it is defined by the approximation of the lines 1-4 extended to R+. Generalizing (2.8) 
to an arbitrary flow, we take 

1 - 1 - 

7= (2 /(n) ) 5 ' ^W'l^l' ^E)) = V (( k -I(u))\e). (2.9) 

If the flow curve is extended by the straight line r = c\ + C27, 7 6 (71, 00) , we obtain 

(/(«), 1^1, e [o,2 7l 2 ], 

J(c2 + ci(i/(u))-3), /(u) G (2 7l 2 ,cx)). (210) 



Here, the coefficients c\ and C2 depend on \E\ and fi(u,E). The viscosity function </? is 
continuous in R+ x [0, 1], if the flow curve is extended in [0, 70] by the dot-dash line, and it 
has the form 

¥>(/(«), \E\ , »(u, E)) = b MlI^)} + |£| , ^ £)) (2.11) 

7(u) 2 

with 6(|.E|, /i(u, -E 1 )) = 2^5 ro, if the flow curve is extended by the dash line in [0,70], Y> 
being a function continuous in R^_ x [0, 1]. 

We note that if ip(I(u), \E\, jx(u, E)) = b\{\E\,n(u,E)), then (2.11) is the viscosity func- 
tion of an extended Bingham electrorheological fluid. 

In the case that the flow curve is extended in [0, 70] by the dot-dash line, the viscosity 
function can be written as follows: 



<p(I(u),\E\,fi(u,E)) =b(\E\, f i(u,E))(\ + I(u))- 1 2 +i>(I(u),\E\,n(u,E)), 



(2.12) 



where A is a small positive parameter. Obviously, for A = the function ip defined by (2.12) 
is the same as the one defined by (2.11). Moreover, for I{u) = we have 

ip(0,\E\,/j,(u,E)) = 00 for (2.11), 

f(fi,\E\,fi(u,E)) = b(\E\,[i(u,E))\~5 +ip(0,\E\,fj,(u,E)) for (2.12). 

Flow problems for fluids with a constitutive equation (2.11) reduce to the solution of vari- 
ational inequalities. Such problems are considerably more complicated than problems for 
fluids with finite viscosity, in particular, for fluids with a constitutive equation as given by 
(2.12). From a physical point of view, (2.12) with a finite, but possibly large viscosity for 
I(u) = seems to be more reasonable than (2.11). 

We will study problems in the case that ip is a continuous bounded function of its argu- 
ments and in the case that <p is singular and the singular part of the function <p> is equal 
to b(\E\, fi(u, E))I(u)~2 . In the first case we assume that (p satisfies one of the following 
conditions (CI), (C2), (C3): 

(CI): tp : (1/1,1/2,2/3) — ► <f{yi, 2/2, 2/3) is a function continuous in R+ x [0, 1], and for an 
arbitrarily fixed (3/2,2/3) G R+ x [0,1] the function (p(.,y 2 ,y 3 ) ■ V\ -»■ ^(2/1,2/2,2/3) is 
continuously differentiable in R + , and the following inequalities hold: 



02 > <p(yi ,2/2,2/3) > a-i (2-13) 
.2/3) 



^(2/1,2/2,^3) + 2 J~ (2/1,2/2,2/3)2/1 > 03 (2.14) 



^ i (2/1,2/2,2/3) 2/1 < 04, (2.15) 

where a^, 1 < z < 4, are positive numbers. 

(C2): (/? : (2/1,2/2,2/3) — * ^(2/1,2/2,2/3) i s a function continuous in R+ x [0, 1], and for an 
arbitrarily fixed (2/2,2/3) 6 ^+ x [0, 1], (2.13) and the following inequality hold: 

[^1,2/2,2/3)^1 - ^(z!) y2,2/3)22](zi - 22) > 03(^1 - 2 2 ) 2 V(zi,z 2 ) G R+. 

(2.16) 

(C3): 99 : (2/1,2/2,2/3) — * V 3 (2/i,2/2,2/3) is a function continuous in R^ x [0,1], and 
for an arbitrarily fixed (2/2,2/3) G ^+ x [0,1], (2.13) holds and the function z — > 
99 (z 2 , 2/2, 2/3)2 is strictly increasing in R + , i.e., the conditions 21,22 £ R+, 21 > z 2 
imply (^(^,2/2,2/3)21 > ^(21,2/2,2/3)22. 

Let us dwell on the physical sense of these inequalities. (2.13) indicates that the viscosity 
is bounded from below and from above by positive constants. The inequality (2.14) implies 
that for fixed values of \E\ and n(u, E) the derivative of the function I(v) — > G(v) is positive, 
where G{v) is the second invariant of the stress deviator 

G{v) = %>{I[v),\E\, fi(u,E))] 2 I(v). 

This means that in case of simple shear flow the shear stress increases with increasing shear 
rate. (2.15) is a restriction on for large values of y\. These inequalities are natural from 
a physical point of view. 

The assumptions (C2) and (C3) indicate that in case of simple shear flow, the shear stress 
must increase with increasing shear rate. 

Moreover, let X(z) = tp(z 2 , j/2, 2/3)2. Assume that the function A is continuously differ- 
entiable in R + . It follows from (C2) that jfe > a 3 for all z G R + . Calculating jfe(z), we 
obtain (2.14) from (2.16). On the contrary, (2.16) follows from (2.14). The assumption (C3) 
indicates that ^(2) > for all z G R + or, equivalently 

^(2/1,2/2,2/3) + 2— (2/1,2/2,2/3)2/1 > 0, yi € R+. 

oyi 

The following assumption is concerned with the function (coefficient) b in (2.11), (2.12): 
(C4): ^ : 2/1,2/2 ^ £>(2/i,2/2) is a function continuous in R + x [0, 1] and in addition 

< 6(2/1,2/2) <a 5 , (2/1, 2/2) G R+ x [0, 1], (2.17) 
05 being a positive number. 
Generally, the continuous function ip is expressible in the form 

m 

( p(I(u),\E\,^u,E))=^e i (\E\,^u,E))p i (I(u)). (2.18) 

i=i 

Polynomials or splines can be used to represent the functions The flow curves obtained 
for various values of E can be approximated with an arbitrary accuracy. We note that (2.18) 
may also be used for an identification of the function ip. 
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In the general case, the coefficients as well as the viscosity function ip depend on the 
temperature and these coefficients can be determined by an approximation of the corre- 
sponding flow curves. 

2.3. General problems. Fig. 3 below gives an example of an electrorheological fluid flow. 



•5*2 



~7 

S 1 = S\S 2 



u 



AU(t) 



Figure 3 



Here, the domain Q of fluid flow consists of three parts CI2, ^3- A fluid flows from the 
part across ^2 m the part f^. Electrodes are placed on parts Tq and Ti of the boundary 
of S72, and an electric field E is generated by applying voltages AU(t) to electrodes at time 
t. Generally it may be a pairs of electrodes and voltages AUi(t) are applied to i-th pair of 
electrodes, i = 1, . . . , k. The boundary S of the domain f2 consists of two parts S± and S 2 . 
Surface forces F = (F±, . . . , F n ) act on 52, and the distribution of velocities u = (ui, . . . , u n ) 
is given on S±. 

The equations of motion and the incompressibility condition read as follows: 



= Ki in Q = flx(0,T), i = l,...,n, 

n 

div u = JT~ = m Q- 



(2.19) 
(2.20) 



Here, Ki are the components of the volume force vector K, p is the density, T a positive 
constant. In (2.19) and below Einstein's convention on summation over repeated index is 
applied. 

We assume that O, is a bounded domain in M. n , n = 2 or 3. Suppose that Si and S2 
are open subsets of S such that S = Si U S2 and Si n S2 = 0. The boundary and initial 
conditions are the following: 



I Si x (0,T) 



u, 



E 



u(.,0) 



u 



in fi. 



Here, Fi and i/j are the components of the vector of surface force F 
unit outward normal v = (ui, . . . , zv n ) to S, respectively. 



1, . . . ,n 



(Fi,. 



(2.21) 
(2.22) 
(2.23) 

, F n ) and the 



We consider Maxwells equations in the following form (see e.g. [8]): 

ldB 

cuv\E + -— = 0, div5 = 0, 
c at 

curlF - = 0, divD = 0. 

c at 



Here £J is the electric field, B the magnetic induction, D the electric displacement, H the 
magnetic field, c the speed of light. One can assume that 



D = eE, B = fiH, 

where e is the dielectric permittivity, \x the magnetic permeability. 

Since electrorheological fluids are dielectrics the magnetic field H can be neglected. 
(2.24), (2.25) give the following relations 

curl£ = 0, 
div(e£) = 0. 

It follows from (2.26) that there exists a function of potential 6 such that 

E = — grad#, 

and (2.27) implies 

div(egrad6>) = 0. 

The boundary conditions are the following: 

6 = Ui{t) on T t , i = l,...,k, 
6 = on r i0 , 

k 

v ■€ grad0 = on S\(\J (r,ur i0 )). 

i=i 



Here Tj and Tjo are the surfaces of the i-th control and null electrodes respectively, and it 
is supposed that Tj, I^o are open subset of S. We assume 



e G ^oo(^), ei < e < e2 a.e. in Q, 
ei, e2 being positive constants. Suppose also that 

Ui(t) € H&iTi), t€[0,T\, i = l,...,k, 

where 

H^Ti), {^ = u| r ., veH 1 ^), v\ sxT - = 0}. 

i 

The space i? 2 (rj) is normed by 

i m = [n £{\\ v \\m(n), veH 1 ^), v\ r _=ip, v\ xr - = 0}. 



Let be a function such that 

9eH\n), e\ r =Ui(t), oL =o, i = i,...,k. 



2.24) 



2.25) 

Then 

2.26) 
2.27) 

2.28) 

2.29) 

2.30) 
2.31) 

2.32) 



2.33) 
2.34) 
2.35) 

2.36) 

2.37) 
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Define a space V and a bilinear form a : H 1 ^) x V — > R as follows: 

F = -{>; G tf 1 ^), ^ =0}, (2.38) 

u^iCTiurio) 

a(v,/i) = / e^^dx, veH 1 ^), heV. (2.39) 
J n dxi dxi 

Consider the problem: find u satisfying 

u G V 

a(u,h) = -a(9,h) heV. (2.40) 

By use of Green's formula it can be seen that, if u is a solution of problem (2.40), then the 
function 9 = u+9 is a solution of (2.29)— (2.32) in the sense of distributions. On the contrary, 
if 9 is a solution of problem (2.29)-(2.32), then u = 9 — 6 is a solution of the problem (2.40). 
Therefore, the function 9 = u + 9 is a generalized solution of problem (2.29)-(2.32). 

Theorem 2.1. Suppose the conditions (2.33), (2.34) are satisfied. Then there exists a 
unique generalized solution 9 of problem (2.29) -(2.32), and the function 9 is represented in 
the form 9 = 9 + u, where 9 satisfies (2.37) and u is the solution of problem (2.40). 

Proof. By virtue of (2.33) the bilinear form a is continuous and coercive in V. Therefore 
there exists a unique solution u of problem (2.40), and the function 9 = 9 + u is a generalized 
solution of problem (2.29)-(2.32). 

Let 9\ = 9\ + u\ and 9 2 = 9\ + u 2 be two generalized solutions of problem (2.29)-(2.32), 
where 9\ and 9 2 satisfy (2.37), and u±, u 2 G V. Then 9\ — 9 2 = w G V, and (2.40) implies 
a(w, w) = 0. Therefore w = 0. ■ 

The functions of volume force K and surface force F in (2.19) and (2.22) are represented 
in the form 

K = K + K e , F = F + F e , (2.41) 

where K and F are the main volume and surface forces, K e and F e are volume and surface 
forces generated by the vector of electric field E. Considering electrorheological fluid as a 
liquid dielectric we present the stress tensor a e = {c e jfc}" fc=1 induced by electric field as 
follows (see [8]): 

a eik = - l -^-(e-p—)S ik + —EiE k . (2.42) 

Taking into account (2.26), we obtain the following formula for the vector of volume force 

K e = (Kel, ■ ■ ■ , Ken), 



K ei 



da eij \E\ 2 



de 1 3 a,, &\ .„ 



dxj 8tt dxi 8ir dxi V dp/ 

The vector of surface forces is given by 

F e = (F el , . . . , F en ), F ei = (JeikVk- (2.44) 

Thus the systems (2.19)-(2.23) and (2.29)-(2.32) are separated, so one can first solve quasi- 
static system (2.29)-(2.32) and then solve the problem (2.19)-(2.23), (2.41). 
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3. Auxiliary results. 

Let Q, be a bounded domain in M. n with a Lipschitz continuous boundary S, n = 2 or 3. 
Let Si be an open non-empty subset of S. We consider the following spaces: 

X = {u\u G H l (tt) n , u\ Sl = 0}, (3.1) 
V = {u\u£X, divn = 0}. (3.2) 

By means of Korn's inequality, the expression 

= (^j I(u)dx) 2 (3.3) 

defines a norm on X and V being equivalent to the norm of if 1 (SI)™. 

Everywhere below we use the following notations: If Y is a normed space, we denote by 
Y* the dual of Y, and by (/, h) the duality between Y* and Y, where / G Y*, h G Y. In 
particular, if / G £ 2 (^) or / G L 2 (0) n , then (/, /i) is the scalar product in L 2 (f2) or in 
L 2 (f2) n , respectively. The sign — ^ denotes weak convergence in a Banach space. 

We further consider three functions v,vi,V2 such that 

veH 1 ^)", Vl eL 2 {n), vi(x) > a. e. in Q, 

v 2 G Loo(n), f 2 (x) G [0, 1] a.e. in 0. (3.4) 

We set v = (v, v\,v 2 ) and define an operator L v : X — > X* as follows: 

(L v (u),h) = 2 / ip(I(u + v), «i, V2)sij{u + v)eij{h)dx u,h£X. (3.5) 

Lemma 3.1. Suppose the conditions (CI) and (3.4) are satisfied. Then the following in- 
equalities hold 

(L v (u) — L v (w),u — w) > fii\\u — w\\x, u,w £ X, (3.6) 

H-Mm) - L v (w)\\ x * < (J.2\\u - w\\x, u,w£X, (3.7) 

where 

li\ = min(2ai, 2a 3 ), /it 2 = 2a 2 + 4a 4 . (3.8) 
Proof. Let u, u> be arbitrarily fixed functions in X and 

h = u — w. (3-9) 

We introduce the function 7 as follows: 

7(i) = / ip(I(v + w + th),v 1 ,V2)e i j(v + w + th)£ij(e)dx, t e [0,1], e G X. 

Jn (3.10) 

It is obvious that 

7(1) - 7(0) = \{L v {u) - L v (w),e). (3.11) 
By classical analysis it follows that 7 is differentiable at any point t G (0, 1). Therefore 

7(1) = 7(0) + ^(0, £€(0,1), (3.12) 
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where 

~dt 

+2-— (I(v + w + £h),v 1 ,v 2 )£km(v + w + ih)e km (h)£ij(v + w + ^h)e ij {e)]dx. 
Cyi 



(0 = / [<p{I(v + w + £h),v 1 ,v 2 )e ij (h)si j (e) 
Jn 

+ £h)eij{e)}dx. 

(3.13) 



Taking note of the inequality 

\£ij(v + w + £h)eij{h)\ <I{v + w + £h)h(h)^ 

and (2.13), (2.15) we get (3.7) as a direct consequence of (3.9)-(3.13). 
Define the function g as follows: 

( \ i ^( a ^Mx),V2(x)), if %fc(a,vi(x),V2(x)) < 0, 
g(a,x) = < ayi % 

[0, if Q^{a,v 1 {x),V2{x)) > 0, 

where a € R+, i6!l. 

Then, taking e = h in (3.13) and applying (2.14) we get 



= jW{v + W + Zh), Vl ,V2)I(h) 



+2-^-(J(fl + w + £h),v 1 ,v 2 )(e ij (v + w + £h)eij(h)) 2 }dx > min(oi, a 3 )||/i|& 

^ (3.14) 

and (3.6) follows from (3.14). ■ 

Lemma 3.2. Let the function ip satisfy condition (C2) and assume that (3.4) holds true. 
Then the operator L v is a continuous mapping from X into X* and 

(L v (u) — L v (w),u — w ) > 2a3 / [I(u + v) 2 — I(w + v) ^] 2 dx, u,w € X. 

Ja (3.15) 

Moreover, 

(L v (u) — L v (w), u — w) = u = w. (3.16) 

Proof. We set u 1 = u + v, w l = w + v. Taking into account that [eij(u l )Eij{w 1 )] < 
I(u 1 )2l(w 1 )^, we obtain 

(L v (u) — L v (w),u — w) = (L v (u) — L v (w),u 1 — w 1 ) 

= 2 / [<p(I(u 1 ),v 1 ,v 2 )I{u 1 ) + (p(I(w 1 ),v 1 ,v 2 )I{w 1 ) 

-(p(I(u 1 ),v 1 ,v 2 )£ij(u 1 )e ij (w 1 ) - (p(I(w 1 ),v 1 ,v 2 )e ij (w 1 )£ ij (u 1 )]dx 

> 2 f [(^(/(u 1 ), ^^/(u 1 ^ - ^(/(ii; 1 ),^,^)/^ 1 )!]^ 1 )! - /(u- 1 )^ 

Jn (3.17) 

Observing (2.16), (3.17) yields (3.15). Now, assume 

(L v (u)- L v (w),u-w) =0. (3.18) 

Then, by (3.15) we have 

I(u + v) = I(w + v) a.e. in 12, ip(I(u 1 ),vi,v 2 ) = (p(I(w 1 ),vi,v 2 ) a.e. in 12. 

(3.19) 
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Taking (2.13), (3.5), (3. 18), (3. 19) into account we get \\u — w\\x = 0. 

The continuity of the mapping L v follows from the continuity of the Nemytskii operator 
(see [24], and also [13], Lemma 8.2 Chapter 2). ■ 

Lemma 3.3. Assume that (C3) is satisfied and (3.4) holds true. Then 

[L v (u) - L v (w),u-w) > 0, u, w £ X. (3.20) 

Moreover, (3.16) is valid, and the operator L v is a continuous mapping from X into X* . 

Proof. Indeed, (3.20) follows from (C3) and (3.17). Assume that (3.18) is valid. Then 
(C3) and (3.17) imply (3.19), and by (2.13) we obtain u = w. ■ 
Define the set U as follows 

U = {h € Loo(n), < h(x) < a 5 a.e. in fi}, (3.21) 

and let v € H 1 (Q) n . For a given constant A > define an operator L x : U x X — > X* as 
follows: 

(L\(h,u),w) = / /i(A + /({; + n))^^ei i (t; + n)ei i (u;)(ia;. (3.22) 
Jn 

Lemma 3.4. For an arbitrary A > and an arbitrarily fixed h £ U the following inequalities 
hold: 

(L x (h,ui) - L x (h,u 2 ), ui - u 2 ) > 0, u 1 ,u 2 £X, (3.23) 
\\L x (h, m) - L x (h, u 2 )\\ x * < a\\ui - u 2 \\x, a = 2a 5 X~^, (3.24) 

\\L x (h,u)\\ x * < an, a x = (^j h 2 dx^j 2 . (3.25) 

Moreover, the conditions 

{h m } C ?7, /i m — > /i a.e. in 

• v • n ^mi 0«i ... 

u m — > -a ml, u m — > « a.e. zn SZ, — ► — — a.e. ?n SZ, z, j = 1, . . . , n, 

dxj dxj (3.26) 

imply 

L x (h m ,u m ) -> L x (h,u) in X*. (3.27) 

Proof. The viscosity function associated with the operator L x (h,.) : u — > L x (h,u) has 
the form 

v(y) = ^(A + j/)-3, (3.28) 
where y plays the role of the second invariant of the rate of strain tensor. We have 
tp(y) + 2^(y)y = h(\ + y)~i[l - (A + y^y] > 0, y G R+, A > 0. 

ay z ^ 3 29 ^ 

The left-hand side of (3.29) represents the derivative of the function g : z — ► <?(z) = i^(z 2 )z, 
z 2 = y. Therefore, the function g is increasing, and (3.23) follows from the proof of Lemma 
3.3. 

By (3.21), (3.22) we obtain 

i 

\(L x (h,u),w)\ < / h(\ + I(v + u))~^ I(v + u)^ I(w)^ dx < ( / h 2 dx) \\w\\ x , 
Jn \Jn J 
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which readily gives (3.25). Moreover, (3.21) and (3.28) yield 



1 



2 



dtp 1 



<p(y) < ~a 5 \-*, -r(y) y < a 5 \—2, yeR+, (3.30) 



dy Kaj a ~ 4 



and hence, observing Lemma 3.1 we get (3.24). 
Assuming (3.26), it follows that 

\\L x (h m ,u m ) - Lx{h,u)\\x* < \\Lx(h m ,u m ) - Lx{h m ,u)\\ x * 

+ \\Lx(h m ,u)-Lx(h,u)\ x *. (3.31) 

Further, 

(L x {h m ,u m ) - Lx(h m ,u),w) = / h rn {(\ + I(v + u m ))~2e ij (u m - u)£ij(w) 

Jn 



+((A + I(v + u m )) 2 — (A + I(v + u)) *)eij{v + u)£ij(w)]dx, 



whence 



mi U m ) - L x (h m ,u)\\x* < 



h 2 m (\ + I(v + u m )) 1 I(u m - u)dx 

V. 

1 



+ {Jh 2 m [{\ + I{v + Um))-* - (A + I(v + u))'^] 2 I{v + u)dx\ 2 . (3.32) 

Obviously, the first term of the right-hand side in (3.32) tends to zero. By (3.26) we have 
I(v + u m ) — > I(v + u) a.e. in Q, and by the Lebesgue theorem we obtain that the second 
term of the right-hand side in (3.32) tends to zero. The second term of the right-hand side 
in (3.31) also tends to zero. Thus (3.27) is satisfied, and the lemma is proven. ■ 

Lemma 3.5. Let Q be a bounded domain in MP, n = 2 or 3 with a Lipschitz continuous 
boundary S, and let the operator B £ C{X, L2(f2)) be defined as follows: 

Bu = div-u. (3.33) 

Then, the inf-sup condition 

inf sup - } B ^^ > Pi > (3.34) 

fieL 2 {n) veX \\v\\x ||H|i 2 (n) 

holds true. The operator B is an isomorphism from V 1 - onto L2(£l), where V 1 - is the orthog- 
onal complement of V in X, and the operator B* that is adjoint to B, is an isomorphism 
from L,2(£l) onto the polar set 

V° = {f ex*,(f,u) = 0, ueV}. (3.35) 

Moreover, 

ll 5_1 ||£(L 2 (n),yi) < (3-36) 

UBT'WcivoMm < ^- (3-37) 

For a proof see in [1]. Lemma 3.5 is a generalization of the inf-sup condition in case that 
the operator div acts in the subspace Hq(Q) (see [6]). This result was first established in 
an equivalent form by Ladyzhenskaya and Solonnikov in [7]. 
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Let {X m }™ =1 , {N m }^ =1 be sequences of finite-dimensional subspaces in X and 1/2(0), 
respectively, such that 

lim inf \\u — z\\x = 0, (3.38) 

lim inf ||w - y|| L2 (m = 0, w G L 2 (fi). (3-39) 

Define the operators B m G £(X m , N^) as follows: 

{B m u,n)= I /idivudx, ueX m , ft G N m , (3.40) 
./n 

and let B^ G C(N m , X^) be the adjoint operator of -B m with (B m u, fx) = {u,B* m [i) for all 
■u G X m and all /x G N m . 

We introduce the spaces V m and by 

F m = {u G X m , (B m tt, n) = 0, /i G iV m }, (3.41) 

^ = {?e^- (?,«) = 0, uGKj. (3.42) 

The following Lemma is valid (see [1]). 

Lemma 3.6. Let {X m }^ =1 , {N m }^ =1 be sequences of finite- dimensional subspaces in X 
and £2(0) and assume that the discrete inf — sup condition (LBB condition) 

inf sup f^"",',^ > > 0, m G N (3.43) 
/ieAr m ugXm ||«||x||m||l 2 (q) 

/ioWs irae. Then the operator B* m is an isomorphism from N m onto V^, and the operator 
B m is an isomorphism from onto N^, where is an orthogonal complement ofV m in 
X m . Moreover, 

UB^Wc^nJ < ^, ll^ll^yx) < i mGN. (3.44) 

Consider a functional * : U x X — > M + of the form 

V(h,u) = f hl(u)*dx h£U, u G X, (3.45) 
Jn 

where U is as in (3.21). 

Lemma 3.7. For an arbitrarily fixed h G U the functional ^(h, .) : u — > ^(h,u) is continu- 
ous in X and the conditions 

{h m } C U, h m h a.e. in $7, u m ^ u in X (3.46) 

imply 

lim inf ^(h m ,u m ) > V(h,u). (3.47) 

Here and below the sign — designates the weak convergence. 
Proof of the Lemma 3.7. Let u rn — ► u in X. We have 

1 . p .1 



hL(u m — u)i dx < J h 2 dx^j 2 ^ y /(u m — u)dx^j 2 . 

lim / /i/(it m - u)3 cfa = 0, (3.48) 



Therefore, 
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and 



/ hl(u m — it) 2 dx > / hl(u m ) 2 dx — hl(u) 2 dx 
Jn Jn Jn 



Consequently, lim ^(h, u m ) = ^>(h,u). 
Let a € [0, 1], u,v € X. Then 



/(to + (1 — a)w) = I(au) + 2a(l — a) ^ £ij(u)eij(v) 
+/((! - a)u) < [al(u)s + (1 - a)/(?;)5] 2 . 



(3.49) 



(3.50) 



Therefore, 

f(/i,att + (l-a)!)) = / W(cm + (1 - a)v)^dx < a^(h,u) + (1 - a)*(/t,' 
which shows that 

^(h, .) : u — > ^>(h,u) is a convex functional in X. (3.51) 
Let now (3.46) be fulfilled. We have 

*(/i m ,u m )= [hl(u m )v + (h m -h)I{u m ) 1 2]dx, (3.52) 
Jn 

(h m - h)I(u m )2 dx < \\h m ~ h\\L 2 (n)\\um\\x- (3.53) 



In view of (3.46), the right-hand side of (3.53) tends to zero as m — > oo. Hence (3.51), 
(3.52), and the continuity of the functional .) imply 



liminf ^f(h m , u m ) = liminf ^(h, u m ) > ^f(h,u), 

and the lemma is proved. ■ 
Remark 3.1. Assume that 

hi < h(x) < ti2 a.e. in Q, 
where hi, hi are positive constants. Then the expression 

hl{u) 2 dx = \\u\\h 



(3.54) 



(3.55) 



(3.56) 



defines a norm on X. Note that this norm is not equivalent to the norm of the space 
Wl(Q) n . However, 



\u\\ p = hl(u)2 dx^j 



(3.57) 



is a norm on X, which is equivalent to the norm of Wp{Q) n for p > 1 (cf., e.g., [15]). 
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4. The stationary problem. 

We consider stationary flow problems of electrorheological fluids under the Stokes ap- 
proximation, i.e., we ignore inertial forces. Such an approach is reasonable, because the 
viscosities of electrorheological fluids are large, and the inertial terms have a small impact. 
We deal with the following problem: find a pair of functions u, p satisfying 

dv d 

^-2 — [i P (I(u),\El»(u,E))e ij (u)}=K i in 17, i = l,...,n, (4.1) 

div u = in n, (4.2) 



u 



Si 



(4.3) 



[-pSij +2<p{I{u),\E\,ii(u,E))e ij {u)]v j = F h i = l,...,n. (4.4) 

S'2 

We assume that 

uem{Si) n . (4.5) 

Then there exists a function u such that 

ueH l (Q) n , u = u, divu = 0. (4.6) 
Si 

Suppose also 

K = (K 1 ,...,K n )eL 2 (Q) n , F = (F 1 ,...,F n )eL 2 (S 2 ) n . (4.7) 
In line with (2.11), we choose the viscosity function <p of the following form: 

¥>(!(«), \E\,»(u, E)) = b MlJ^)} + mu) , lElKU: E )\ (4.8) 

7(u)2 

where ip is a function satisfying one out of the conditions (CI), (C2), (C3) with <p replaced 
by ip, and b satisfies (C4). We refer to the fluid with the viscosity function ip defined by 
(4.8) as a generalized Bingham electrorheological fluid. 

Define a functional J on the set X x X and an operator L : X — > X* as follows: 

J(v,h)=2 I b(\E\,/j,(u + v,E))I(u + h)^dx, v,heX. (4.9) 



(L(v),h)=2 / ip(I(u + v), |^|,/x(tf + v,E))eij(u + v)£ij(h)dx, v,heX. 
Jn 



(4.10) 



We use the notations 



(K,h)= [ Kihidx, (F,h)= [ Fihids, heX. (4.11) 
The following assertion holds. 

Theorem 4.1. Lei £7 &e a bounded domain in M. n , n = 2 or 3, mi/t a Lipschitz continuous 
boundary S. Assume (4.6), (4.7) are satisfied and (u,p) with u = u + v is a regular solution 
of (4.1) -(4.4), where the viscosity function ip is defined by (4.8) with ip meeting one out of 
the conditions (CI), (C2), (C3) (ip replaced by ip) and b satisfying (C4). Then 

v e V, (4.12) 

J(v,h)-J(v,v) + (L(v),h-v)>(K + F,h-v), heV. (4.13) 
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Proof. Let h = (hi,..., h n ) G V. We multiply (4.1) with hi — vi, sum over i and integrate 
over Q. By Green's formula and (4.4), (4.8), we obtain 

2 / b(\E\,^(u + v,E))I(u + v)~*£ij(u + v)eij(h — v)dx 
Jn 

+(L(v),h-v) = (K + F,h- v), heV. (4.14) 

We use the relations 

£ij(h - v) = Eij(u + h) — £ij(u + v), £ij(u + v)£ij(u + h) < I(u + v)*I(u + 

(4.15) 

so that the first addendum of the left-hand side of (4.14) is majorized by J(v, h) — J(v,v). 
Then, (4.14) implies (4.13), and the theorem is proved. ■ 
Let v be a solution of the problem (4.12), (4.13) such that 

I(u + v) / a.e. in n. (4.16) 

We replace h in (4.13) by v + Xh, A > 0, from which 

A _1 [J(v,v + A/t) - J(v,v)] + {L(v),h) > (K + F, h). (4.17) 

For A — > we get 

dJ 



. )h (v,v),h) +(L(v),h)>(K + F,h), heV, [IAS) 

where j^(v,v) is the partial Gateaux derivative of the functional J with respect to the 
second argument 



— (v,v),h \ = 2 J b(\E\,fi(u + v,E))I(u + vyhi j (u + v)e ij (h)dx. 



Replacing h by — h in (4.18) we obtain 

dJ 



(4.19) 



;)h (v,v),h) +(L(v),h) = (K + F,h), heV, (4.20) 



and Lemma 3.5 gives 

dJ 



;)h (v,v)+L(v)-K-F = B*p, P eL 2 (n), (4.21) 



that is 

dJ 



;)h (v,v),hj +(L(v),h)-(B*p,h) = (K + F,h) VheX. (4.22) 

It follows from (4.22) that the pair (u,p) with u = u + v is a solution of (4.1)-(4.4) in the 
sense of distributions. Thus, we have proved the following: 

Remark 4.1 If v is a solution of (4.12), (4.13) that satisfies (4.16), then there exists a 
function p G £2^) such that the pair (u,p) with u = u + v is a solution of (4.1)-(4.4) in 
the distributional sense. In view of this and Theorem 4.1 it is reasonable to refer to the 
function u = u + v as a generalized solution of (4.1)-(4.4). 
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5. Problem for the fluid with constitutive equation (2.12). 
5.1. Existence theorem. We define the following functional on X x X: 

J x {v,h) =2 f b(\E\,ii(u + v,E))(\ + I(u + h))*dx, A > 0. (5.1) 
Jn 

Obviously, J\(v, h) = J(v, h) for A = 0. Note that the functional J\ is Gateaux differentiable 
in X with respect to the second argument for A > 0, but not for A = 0. 
The partial Gateaux derivative ^ is given by 

(^-(v,h),wj = 2 J b(\E\, l i(u + v,E))(X + I(u + h))~h ij (u + h)£ ij (w)dx, 

v,h,w£X, A>0. (5.2) 
Consider the following problem: find v\ such that 

vx G V, (5.3) 
(-Q±(v x ,v x ),w} +(L(v x ),w) = (K + F,w), w G V. (5.4) 

Lemma 3.5 implies that if v\ is a solution of (5.3), (5.4), then there exists a function px 
such that the pair (vx, Px) is a solution of the following problem: 

(vx,Px) G X x L 2 (Q), (5.5) 
{^(vx,v x ),w) +(L(vx),w)-(B*px,w) = (K + F,w), w G X, (5.6) 
(Bv x ,q)=0, qeL 2 (Q). (5.7) 

We remark that (5.5)— (5.7) represent the flow of an electrorheological fluid with the consti- 
tutive equation (2.12). We seek an approximate solution of the problem (5.5)-(5.7) of the 
form 

( v m,Pm) G X m X N m , (5.8) 

d J 

(-^(fm^m),™) + (L(v m ),w) - (B^p m ,w) = (K + F,w), w G X m , (5.9) 

{B m v m ,q)=0, qeN m , (5.10) 

where X m and N m are finite dimensional subspaces in X and L 2 , (fi), respectively, and B m 
is defined as in (3.40). 

Theorem 5.1. Suppose that the conditions (C4), (4.6), (4.7) are satisfied and the function 
tp meets one of the conditions (CI), (C2), (C3) (ip replaced by ip). Let {X m }, {N m } be 
sequences of finite- dimensional subspaces in X and L 2 {p), respectively, such that (3.38), 
(3.39), (3.43) hold and 

X m C X m+1 , N m C N m+1 , m G N. (5.11) 

Then, for an arbitrary A > there exists a solution [yx, Px) of (5.5) -(5.7). Moreover, 
for m G N and A > there exists a solution (v m , p m ) of (5.8)-(5.10), and a subsequence 
{{vk,Pk)} ca - n be extracted from the sequence {(v m ,p m )} such that — 1 vx in X , p^ — px 
in L 2 (£l). 
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Proof. It follows from (3.41), (5.8)-(5.10) that solution of the problem 

v m eV m , (-7^-(v m ,v m ),u?J + (L(v m ),w) = (K + F,w), w G V m . 

Taking (2.17) into account we obtain 

(^.),.)|-2| / m ,Kn + e,E)) S f + ^h x 
K d h ' Jn (A + /(n + e))5 

< 2 / b(\E\,n(u + e, E))I(e) 2 dx < ci\\e\\ x , A > 0, eel, 
Jn 



(5.12) 



where 

ci = 2a5(mes Q)? . 
By (2.13), (4.6), (4.7), (4.10), and (5.13), for an arbitrary e G X we get 
'0J A , 



(5.13) 
(5.14) 



*(e) 



<9/i 



e,e),e) +(L(e),e) - (if + F,e) > 2ai||e||^ -c||e||x, e G X, A > 0, 



(5.15) 



c 

2ai - 



giving z(e) > for \\e\\x > r 

From the corollary of Brouwer's fixed point theorem (cf.[5]) it follows that there exists a 
solution of (5.12) with 

\\v m \\x < r, \\L(v m )\\ x * < c 2 , m G N, (5.16) 

where the second inequality follows from (2.13) and (4.6). For an arbitrary / G X* we 
denote by G/ the restriction of / to X rn . Then G/ G X^, and by (3.42), (5.12) we obtain 

G(^(v m ,v m ) + L(v m )-K-F) eV°. (5.17) 
Therefore, there exists a unique p m G N m (see Lemma 3.6) such that 



B* mPm = G(-^{v m , v m ) + L(v m ) -K-f). 



(5.18) 



Thus the pair (v m , p m ) is a solution of (5.8)-(5.10). Due to (2.17), (4.7), (5.16) and Lemmas 
3.4, 3.6 we get 



\\Pm\\L 2 (n) < c, m G N. 
By (5.16), (5.19) we can extract a subsequence {v^^p^} such that 

v v vq in X, 
v v — > t>o in L2(^) ra and a.e. in f2, 
p v ^Po inL 2 (n), 
^K)^X m A*, 



3/i 



K,Ur,)^Xi in X*. 



(5.19) 

(5.20) 
(5.21) 
(5.22) 
(5.23) 

(5.24) 



Let 770 be a fixed positive number and w G X Vo , q G A^ . Observing (5.20), (5.22)-(5.24) 
we pass to the limit in (5.9), (5.10) with m replaced by 77, and obtain 

(Xi+X- B*p ,w) = (K + F,w), w£X V0 , 



J 

Jn 



q div vq dx = 0, q G iV„ . 



(5.25) 
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Since rjo is an arbitrary positive integer, by (3.38), (3.39) 

Xi + X ~ B* Po = K + F, 
div^o = 0. 

We present the operator L(v) in the form 

L(v) = L(v,v), 



(5.26) 
(5.27) 

(5.28) 



where the operator (v, w) — ► L(v, w) is considered as a mapping of X x X into X* according 
to 



(L(v,w),h) = 2 / ^(I(u + w),\E\,fi(u + v,E))£ij(u + w)eij(h)dx. (5.29) 
Jn 



We get 



X v (w) = (-^( v v> v v) + l ( v t),v v ) ~ ~^( V V' W ) ~ L(v v ,w),v v -w"j, w eX. 
Lemmas 3.3, 3.4 imply 



We have 



X v (w) > 0, i)£N, w £ x. 



dJ\ dJ x , . 



(5.30) 
(5.31) 



<2[ / [b(\E\, n(u + v v ,E))- b{\E\, [i{u + v , E))] 2 dx] 2 . (5.32) 



(2.17), (5.21), (5.32) and the Lebesgue theorem give 



-(*;,,, w) -> — (t;o,?«) m A . 



Likewise we obtain 



L(v„,w) — > L(vq,w) in X* 



Taking into account that (B v v v ,p v ) = 0, by (5.9), (5.20), (5.22) we obtain 
d J 

(-q^(v v ,v v ) + L(v v ),Vr^ = {K + F,v v ) -► (K + F,^), 

and 



(5.33) 
(5.34) 

(5.35) 

(5.36) 



lim (-^( v V' v v) + L ( v r ] )i w ) ~ ( B *Po,w) = (K + F, w), w G X. 
Observing (5.33)-(5.36) and passing to the limit in (5.30), by (5.27), (5.31) we get 

K + F - -7^(v ,w) - L(vq, w) + B*p , v - w J > 0, w € X. (5.37) 



We choose w = vq — 7/1, 7 > 0, h € X, and consider 7^0. Then, Lemmas 3.3, 3.4 give 



dJ 

K + F- -q£(vo, vq) - L(v ,v ) + B*p , h)>0. 



(5.38) 



This inequality holds for any h £ X. Therefore, replacing h by —h shows that equality 
holds true. Consequently, the pair (v\, p\) with v\ = vq and p\ = po solves (5.5)-(5.7). The 
theorem is proved. ■ 
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5.2. On the uniqueness of the solution. Let v\, w x be two solutions of (5.5)-(5.7) and 

e = w x - v x . (5.39) 

Define a function r] as follows: 

ri(t) = (-^(vx + te, v x + te),e) + (L(« A + te),e), *€[0,1]. (5.40) 
It follows from (5.6) that 

?7(1) - r)(0) = 0. (5.41) 

Assume that the function /j is defined by (2.7). Suppose also that the functions b and tp are 
continuously differentiable and in addition 

j^(yi,\E(x)\,y 3 ) y{ < c, ( yi ,y 3 ) € R+ x [0, 1], xefi. (5.42) 

Note that (5.42) is a restriction on the behaviour of the function at large values of y\. 
Under the above conditions the function r] is differentiable, and we have 

1/(0) = ^(0, e € (0, 1). (5.43) 



Here 



where 



i=i 



71 (0 =2 y ^(|£7|,/ i (« + « A + ee, J E7))/ { (A + /(fi + i; A + ^))-3 

xe„ (u + ^a + ^e)ejj (e)dx. (5.44) 



, _J al +u + v x + t,e h \ r/ e £, \ 

f_ U/n + \u + v A + Ce| ' \E\JW"- Wcty/n + |u + v x + £e| ' 
-(ai/n + |tt + v\ + £e|)~ 2 |w + v x + ^"e| — 1 



/ E \ ~ 

y.{{ui + v xi + £ei)ei)(al + ii + v x + £e, ^) Rn , (5.45) 

72(0=2 f b(\E\,ii(u + v x + Se,E))[-(\ + I(u + v x + Ze))-* 
Jn 

X£km (u + v x + e km (e) Sij(u + v x + £e)Eij(e) + (A + I(u + v x + £e))~2 J( e )] 

(5.46) 

73(0=2 / ^(I(u + v x + te),\E\,fi(u + v x + t;e,E))fz 
Jn oyz 

X£ij(u + v x + £e)£ij(e)dx, (5-47) 

74(0=2 / mi(u + v x + tie),\El f i(u + v x + te,E))I(e) 
Jn 

+2^(I(u + v x + f e), M« + u A + £e, #)) 
X£km(u + v x + ie)e km (e)e ij (u + v x + £e)sij(e)] dx. (5.48) 
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By using (4.15) it is easy to see that 

72(0 > 0, (5.49) 

and Lemma 3.1 implies 

74(0 > Mi||e||x) /"i = min(2ai,2a 3 ). (5.50) 

(5.45) yields 

\M < 4 — 7= — n- (5-51) 

We denote 

db 

b = sup —(\E(x)\,y 2 ) , jft €[0,1], xea (5.52) 
By (5.44), (5.51), and (5.52) we obtain 

|7i(0l < 86o||e||x( / \e\ 4 dx) 4 ( f (a^ + \u + v x + £e|)" 4 cfe) 4 < ci&o||e||^. 

Wn ' Wn ' (5.53) 

Here 

ci = 8cc, (5.54) 

where c is the constant of the inequality 

|||e|||L 4 (fi) < c||e||x, (5.55) 

and 



sup( / {a^i+\u + vx + te\) 4 dx)*, f€(0,l), ||e|| x < 2r, 



(5.56) 



r being the constant of (5.16). 

(5.42), (5.47), (5.51), and (5.53) yield 

< 8c||e|| x ( [ |e| 4 dxY ( [ {a^n~ + \u + v x + Ce|)~ 4 dx) 3 < cyc\\ef x . 

Wq ^ y (5.57) 

Assume that 

Mi - ci(6 + c) = c > 0. (5.58) 

Then > c ||e||^, and by (5.41), (5.43) we obtain that e = 0. Thus, we proved the 

following: 

Theorem 5.2. Suppose that the conditions (CI) (V replaced by ip), (C4), (4.6), (4.7), 
(5.42), (5.58) are satisfied. Then there exists a unique solution of (5.5) -(5.7) in the ball 

d r = {ueX, \\u\\ < r = a^\\K + F\\ x *}- 

Note that in the case that the values of the function \E\ are small the constants 60 an d c 
are small, (5.58) is satisfied and there exists a unique solution of (5.5)-(5.7). 
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6. Variational inequality for the extended Bingham electrorheological 

fluid. 

We now consider a problem on stationary flow of the extended Bingham electrorheological 
fluid. The constitutive equation of this fluid is the following: 

<p(I(u), \E\,n(u, E)) = b (\ E \^{u,E)) + MjBU(n) E) y (6 1} 

We deal with the problem (4.1)-(4.4) and assume that (4.5) and (4.7) are satisfied. Then, 
according to Remark 4.1 the generalized solution of our problem is u = u + v , where u is a 
function satisfying (4.6) and v is a solution of the problem 

v G V, (6.2) 
J(v,h)-J(v,v) + (L 1 (v),h-v)>(K + F,h-v), h G V. (6.3) 

Here, J is the functional given by (4.9) and the operator L\ : X — > X* is defined as follows 

(Ly(v),h) = 2 / b 1 (\E\,ii(u + v,E))e ij (u + v)e ij (h)dx, v,heX. (6.4) 
Jn 

The function b\ is subject to the following condition: 

(C5): b\ : (2/1,2/2) — ► b± (2/1,2/2) is a continuous function on M + x [0, 1] and satisfies 

a 6 < 61(2/1,2/2) < a 7 , (2/1,2/2) et+ x [0,1], (6.5) 

with positive constants a6 and 07. 

We approximate the functional J by J x as given by (5.1). Replacing J by J^, by analogy 
with the reasoning in the proof of Theorem 4.1 we obtain the following problem: find v x 
such that 

v x G V, (6.6) 
{^(v x ,v x ),w)+(L 1 (v x ),w) = (K + F,w), weV, (6.7) 

with ^ given by (5.2). 

Theorem 6.1. Suppose that conditions (4.6), (4.7), (C4), (C5) are satisfied. Then there 
exists a solution v of (6.2), (6.3). Moreover, for an arbitrary A > there exists a solution 
of (6.6), (6.7), and there exists a function p x such that the pair (v x , p x ) is a solution of the 
problem 

(v x , Px ) e X x L 2 (n), (6.8) 

(-^(v x ,v x ),w)+(L 1 (v x ),w)-(B*p Xl w) = (K + F,w), w G X, (6.9) 

(Bv x ,q)=0, qeL 2 (Q). (6.10) 

A subsequence can be extracted from the sequence {v x }, again denoted by {v x }, such that 

v x — v in X and v x — > in L2(fJ) n as A — > 0. (6-H) 

If /(u + v) 7^ almost everywhere in Q, then the functional h — > J(t>, /i) is Gateaux differ- 
entiate at the point v, and there exists a function p G £2^) suc/t i/iai i/ie pair (-y, p) is a 



2d 



solution of the problem 

v £ V, p £ L 2 (n), (6.12) 
(—(v,v),h)+(L 1 (v),h)-(B*p,h) = (K + F,h), heX, (6.13) 

with |^ given by (4.19). 

Proof. The existence of a solution v\ of (6.6), (6.7) follows from Theorem 5.1 as well as 
the existence of a function p\ such that (6.8)-(6.10) hold. It is inferred from the proof of 
Theorem 5.1 (see (5.15)) that 

v\ remains in a bounded set of V independent of A. (6-14) 

Therefore, from sequence {v\} we can select a subsequence, again denoted by {vx}, such 
that 

v \ — v in X as A — > 0, (6.15) 

v \ — ► v in Z^fi)™ and a.e. in fi. (6.16) 

For /i G V we introduce 

Z A = (Li(v x ), h-v x ) + Jx(v x , h) - Jx(v x , v x ) — (K + F, h — v x ). (6.17) 
Using (6.9), we see that 

Zx = - (jQ^(vx,vx),h - v^j + J x (v x ,h) - Jx(v x ,v\)- (6.18) 

It follows from (5.1), (5.2) and Lemma 3.4 (cf. (3.23)) that for an arbitrarily fixed w 6 X 
the functional u — > Jx(w,u) is convex. Therefore 

Z A > 0. (6.19) 

(C5), (6.16) and the Lebesgue theorem give 

6i(|S|,^(« + « A ,£7))e ij (/i)^6i(|£;|,/i(« + u), J E;))e ij -(/i) in L 2 (0). (6.20) 

(6.15), (6.20) imply 

lim(Li(u A ),/i) = (Li(i;),/i). (6.21) 



We have 



where 



(L 1 (v x ),vx) = A 1X + A 2X , (6.22) 



A 1X = 2 f hdEl^iu + vx^^jiu^Mdx, (6.23) 

A 2X = 2 [ h(\E\,Li(u + vx,E))I(vx)dx. (6.24) 

(6.20) still holds true if the function h is replaced by u. Consequently, (6.15) implies 

limA lx = 2 f b^El^u + v^E^eij^Si^v). (6.25) 
Jn 

It follows from (6.5) and (6.16) that 

[h(\E\,p(u + vx,E))]*w^[h(\E\,p(u + v))]*w inL 2 (H), u)ei 2 (Sl), 
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so that (6.15) yields 

lim / [b 1 {\E\^{u + v X: E))]^£ ij {v x )wdx = I [b 1 {\E\,ii{u + v,E))\^£ ij {v)wdx, 
Jn Jn 

w e L 2 {yt). 

Therefore, 

[hdElfiiu + v^E))]^^) - [b^El^u + v^^e^v) in L 2 (Q). 

(6.26) 

(6.24) and (6.26) give 

liminf A 2X > 2 / + v, E))I(v)dx. (6.27) 

By (6.22), (6.25), and (6.27) we obtain 

hminf^K),^) > (Li(v),v). (6.28) 
(5.1), (6.16) and the Lebesgue theorem give 

limJ x (v x ,h) = J(v,h). (6.29) 

Setting 

we have 
where 

and 



b x = b(\E\,fi(u + v x ,E)), b = b(\E\, f i(u + v,E)), (6.30) 
I x = I(u + v x ), I = I(u + v), (6.31) 



J\(v X ,v x ) = J x (v, x ) + B 1X , (6.32) 

= 2 f (b x -b )(X + I x )Ux, (6.33) 
Jn 



D 



1A 



\B 1X \<2(^(\ + I x )d x y (fjb x -b \ 2 d x y. (6.34) 
(6.15), (6.16) and (2.17) imply 

li mJ BiA = 0. (6.35) 

(4.9) and (5.1) yield J x (v,v x ) > J(v,v x ), whence 

liminf J x (v, v x ) > liminf J{v, v x ). (6.36) 

(6.15) and Lemma 3.7 yield 

liminf J(v,v\) > J(v,v). (6.37) 
(6.32), (6.35), (6.36), and (6.37) give 

liminf J x (v x ,v x ) > J(v,v). (6.38) 

(6.10) , (6.15) imply (6.2), and using (6.17), (6.19), (6.21), (6.28), (6.29), (6.38) we obtain 
(6.3). 

It follows from Remark 4.1 that if I{u + v) / almost everywhere in 0, then there exists 
a function p such that (6.12), (6.13) hold. ■ 

Remark 6.1. Assume that in (6.3), (6.7), (6.9) the operator L\ is replaced by the 
operator L defined by (4.10), and the function -ip meets one of the conditions (CI), (C2), 
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(C3) with <p replaced by ip. Then, by Theorem 5.1 for an arbitrary A > there exists a 
solution (v\, p\) of (6.8)-(6.10), and a subsequence {v\} can be extracted satisfying (6.15), 
(6.16). 

However, (6.15), (6.16) do not imply liminf(L(^), v\) > (L(v),v) (compare with ((6.28)), 
and we cannot assert that v is a solution of (4.12), (4.13). In the next section we prove 
the existence of a solution of (4.12), (4.13) under conditions which are more restrictive than 
those of Theorem 5.1. 

7. General variational inequality. 

We assume that the function \i in the operator L defined by (4.10) is replaced by a 
function fj,i such that 

u m ^u in X => yLi(u m ,E) -> fxi(u, E) in L oo (0). (7.1) 

According to (2.7) we may define [i\ as follows: 

_ / aI + Pu{x) + u pi + PE(x) \ 2 
AX) ~\a^i+\Pu{x)+uV (3^i+\PE{x)\)^ [ ' 

where a, (3 are small positive constants, I is a vector with components equal to one, and P 
an operator of regularization given by 

Pu(x) = u{\x - x'\)u(x')dx' , x£~Ti, (7.3) 

where 

uj e C°°(R + ), supp uj = [0, a], u{z) > z£R + , 

uj(\x\)dx = 1, a is a small positive constant. (7.4) 



/ 



In (7.3) we assume that the function u is extended to M n . In case that Pu(x) / a.e. in Q 
we may choose a = 0, if PE{x) / a.e. in Q we may choose (3 = 0. 

For the function \i\ condition (7.1) is satisfied. From the physical point of view (7.2) 
means that the value of the function \i\ and therefore the viscosity of the fluid at a point 
x depends on the angle between the vectors of velocity and electric field strength at points 
belonging to some small vicinity of the point x, implying that the model is not local. 

This seems to be natural, since electrorheological properties of the fluid are linked with 
the presence of small solid particles in the fluid. The mean dimension of these particles can 
be taken as the regularization parameter a. 

In the case under consideration the operator L is defined as follows: 

(L(v),h) =2 / tp(I(u + v),\E\, fii(u + v, E))eij(u + v)eij(h)dx, v,heX. 

Jil ' (7.5) 

We assume also that the following condition of uniform continuity of the function ip holds: 
(C6): for an arbitrary 7 > there exists e > such that the conditions 

2/3.2/3 G [0,1], Wz-t/i\<e, 2/1,2/2 £«+ 

imply 

1^(2/1,2/2,2/3) - ^(2/1,2/2,2/3! < 7- 
The function \x can as well be replaced by the function fii in the functionals J and J\. 
In the following theorem this is not assumed, although it is also valid in this case. 
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Theorem 7.1. Suppose that the conditions (C4), (4.6), (4.7) are satisfied and assume that 
the function ip meets one of the conditions (CI), (C2), (C3) (^ replaced by ip) and that 
(C6) holds. Further assume that the function fii meets condition (7.1), and £/te operator L 
is given by (7.5). Then, for an arbitrary A > there exists a solution (v\, p\) of (5.5)- 
(5.7), and there exists a solution v of (4.12), (4.13). A subsequence can be selected from the 
sequence {v\}, again denoted by {v\}, such that 

v\ — ^ v in X and v\ — ► v in L 2 (S1)™ as A — > 0. 

Proof. 1) The existence of a solution (v\, p\) of (5.5)-(5.7) follows from Theorem 5.1, 
and it is inferred from the proof of this theorem (see (5.15)), that v\ remains in a bounded 
set of V independent of A. Therefore, from the sequence {v\} we can select a subsequence, 
again denoted by {v\}, such that 

v x ^v mV as X^O, (7.6) 

v \ — ► v in L2(H) n and a.e. in fi. (7.7) 

For every v\ we define a functional *a as follows: 

^x(v) = Jx(vx,v) + Mv x ,v)-(K + F,v), v £ V, (7.8) 

where 

${v x ,v) = y (y v(c,i^Ui^+^a,^)k)^- (7.9) 

Consider the problem: find a function v satisfying 

»eV, ^\{v) = min *aW- (7.10) 

If w is a solution of (7.10), then we have 

vGV, -y x (v + th) =(-^(v x ,v),h) + (L(vx,v),h)-(K + F,h)=0, h £ V. 

at t=o \ uv / y\\ 

Here, L(v x , v) is the Gateaux derivative of the functional <&(vx, .):«—> $>(v\, u) at the point 
v, i.e. 

— (v x ,v) = L(v x ,v). 
The operator L(vx, .) : I 3 « -* L(vx, u) G X* has the form 

(L(vx,u),h) = 2 / V'(^(« + ^)J^| ) /^i(« + WA ) ^))eji('W + «)eii(^)dx, u,h£X, 

J n (7.12) 

and (7.5) yields 

L(u,u) = L(v). (7.13) 

It follows from (5.3), (5.4) that the function v = vx is a solution of (7.11). By means 
of Lemmas 3.1-3.4 the functional *S>x is strictly convex. Therefore, there exists a unique 
solution v = vx of (7.11), and problems (7.10) and (7.11) are equivalent. 
(7.10) implies 

Jx(v x , h) + $(v A) h) - J x (vx, vx) ~ Hvx, v x ) >(K + F,h- v x ), heV. 

(7.14) 
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2) It follows from (7.6), (7.7) and the proof of Theorem 6.1 (see (6.29), (6.38)) that 



We have 



where 



lim Jx(v x ,h) = J{v,h), liminf J\{v x , v\) > J{v,v) 

&{v\,h)= / fxdx, §(v,h) = / fdx, 
Jn Jn 



(7.15) 
(7.16) 



rl{u+h)(x) 

h(x) = / ^,\E(x)\,^(u + v x ,E)(x))d^, 

J o 

rl(u+h)(x) 

f(x) = / ^,\E(x)\,^(u + v,E)(x))d^. 

Jo 

(7.7) and (C6) imply f\—*f almost everywhere in $7, and (2.13) yields 

\fx\ < a 2 I(u + h). 
Thus (7.16) and the Lebesgue theorem give 

lim<I>(t; A , h) = $(v, h). 

It is obvious that 

Hv x ,v x ) = &(v,x) +a x , 
Jn L Jo 



(7.17) 
(7.18) 



-^Elu^u + ViE^d^dx. 



(7.1), (7.6) and (C6) imply 



"A < Px 



\ I(u + v x)dx, 
Jn 



and lim/?A = 0. Therefore limaA = 0. 

The functional : u — > $(v,u) is continuous in X. Indeed, let u m —>■ u in X. We 

have 



< a 2 



•I(u+Um) 

~{u+u) 

I (I(u + u m ) - I(u + u))ete 
Jn 



r . rl(u+u m ) 

<Z>(v,u m ) - <S>(v,u) =1(1 ^(£,|£7|,/ii(u + 

Jn v Juu+u) 



dx 



(7.19) 



and the right hand side of this inequality tends to zero asm^oo. 

By Lemmas 3.1, 3.2, 3.3 the functional &(v, .) is convex in X. Therefore, <3?(t>, .) is lower 
semi-continuous for the weak topology on X, and (7.6), (7.18) imply 

liminf &(v x ,v x ) > $(v,v). (7.20) 

By (7.15), (7.17), (7.20) we pass to the limit as A ^ in (7.14). This gives 

J(v,h) + $(v,h)-J(v,v)-$(v,v)>(K + F,h-v), heV. (7.21) 
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Taking into account (7.21) and the convexity of the functional J(v, J(v,u), we get 

J(v, v) + v) — (K + F, v) 
< J{v, (1 - 9)v + 0/i) + (1 - 0)u + 0/t) - (if + F, (1 - 0)u + Oh) 
< (1 - 0) J(v, v) + 0J(v, /i) + (1 - 0)u + 0/t) - (K + F, (1 - 0)u + 0/t), G (0, 1), 
whence 

— — ^— - + J(v, /i) - J(u, u) - (K + F, h - v) > 0. 

(7.22) 

For 9 -> we get (4.13), with the operator L defined by (7.12), (7.13). ■ 

Remark 7.1. We have shown that (7.21) implies (4.13). Let us show that (4.13) yields 

(7.21), that is problems (4.12), (4.13) and (4.12), (7.21) are equivalent. 
Let (4.12), (4.13) be valid. Obviously, 

J(v, h) + h) - J(v, v) - <Z>(v, v) - (K + F,h-v) 

= J(v, h) - J(v, v) + (L(v, v), h - v) + h) - <3?(v, v) 

-(L(v,v),h-v)-(K + F,h-v) (7.23) 

The functional u — > $(v,u) is convex, whence 

/i) - u) - u), /i - u) > 0. (7.24) 

(4.13), (7.13), (7.23), (7.24) give (7.21). 

Remark 7.2. By comparing Theorems 5.1 and 7.1 we observe that a solution of the oper- 
ator equations for a fluid with bounded viscosity function (2.12) exists under less restrictive 
conditions than the conditions for the existence of a solution of the variational inequality 
(4.13) for a fluid with unbounded viscosity function (2.11). In addition, such an important 
characteristic of the flow as the function of pressure is defined for (2.11) only in case that 
(4.16) holds, i.e., when the viscosity function dees not take infinite values. But in this case, 
the variational inequality reduces to operator equations as outlined in Section 4. 

Moreover, from a physical point of view a fluid with finite viscosity (2.12) seems to be 
more reasonable than a fluid with unbounded viscosity (2.11). 

8. Problems with given function fi. 

In the case that the distance between the electrodes is small compared with the lengths 
of the electrodes, one can assume that in between the electrodes the velocity vector is 
orthogonal to the vector of electric field strength, and the electric fields strength is equal to 
zero in the remaining part of the domain under consideration. 

In this case one reckon that fi(u, E) is a known function of x, so that the viscosity functions 
(2.11) and (2.12) take the form 

V (I(u),\E\,x) = ^^ + Ml(y'),\E\,x), (8.1) 

I(U) 2 

ip{I{u),\E\,x) =e(\E\, x)(A + /(«))-3 + \E\, x), (8.2) 

and the constitutive equation is defined by (2.1). 

Dependence of the viscosity function on x in (8.1), (8.2) is connected with the anisotropy 
of the fluid. If the direction of the velocity vector at each point x at which E(x) ^ is 
known, then the viscosity functions (2.11) and (2.12) transform in relations (8.1), (8.2). 

We assume the function ipi to satisfy 
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(CO): for almost all x G ft the function ipi(-, -i x ) '■ (2/1,2/2) — ► 01 (2/1, 2/2, a?) is continu- 
ous in R+, and for an arbitrarily fixed (3/1 , J/2) G R+ the function i/ji(yi, 2/2, •) : £ — > 
■01 (yi> 2/2, is measurable in f2. 

We also suppose that for almost all x G SI and all 3/2 G R+, the function 
V-'iG, 2/2, a?) : — > ipi(yi,y2,x) satisfies one of the following conditions (Cla), (C2a), 
(C3a): 

(Cla): ipi(.,y2,x) is continuously differentiable in R+ and the following inequalities 
hold: 



0-2 > tpi(yi,y2,x) > ai, (8.3) 
.t/2,a 



ipi(yi,y2,x) + 2-Q—(yi,y 2 ,x) > a 3 (8.4) 



(yi,y2,x / 



yi < a A - (8.5) 



(C2a): (8.3) is fulfilled and for an arbitrary (z\,z 2 ) G R+ the following inequality is 
valid: 

[ipi(zf,y 2 ,x)z 1 - ipi(z 2 ,y 2 ,x)z 2 ](z 1 - z 2 ) > a 3 (zi - z 2 f . (8.6) 

(C3a): (8.3) is fulfilled and the function z — > ipi{z 2 , y 2 , x)z is strictly increasing in R + , 
i.e., the conditions 21,22 £ z i > z 2 imply tpi(zf,y 2 ,x)zi > ijji(z 2 l ,y2,x)z 2 . 

(Cla), (C2a), (C3a) are analogies of conditions (CI), (C2), (C3), and an analog of (C4) 
is the following condition: 

(C4a): for almost all x G fl the function e(.,x) : y — ► e(y,x) is continuous in R + and 
for an arbitrarily fixed y G R + , the function e(y, .) : x — > e(y,x) is measurable in f2 
and 

0<e(y,x)<a 5 . (8.7) 
Define functional F and 1a, A > 0, as follows: 

y(«) = 2 /" e(\E\,x)I(u + u)*dx, ueX, (8.8) 

y A («) = 2 / e(|£;|,a:)(A + /(« + u))3da:, u e X. (8.9) 
Define also an operator L 2 : X — ► X* by means of 

(L 2 (u),h) =2 ^ 1 (I(u + u),\E\,x)e i j(u + u)e i j(h)dx, u,heX. (8.10) 



Consider the following two problems: 

Problem 1. Find a pair of functions (v x , p x ) such that 

vaGX, p x eL 2 (Q), (8.11) 

(^K), /») + (L 2 (v x ), h) - (B*p x , h) = {K + F,h), heX, (8.12) 
(Bv x ,q)=0, qeL 2 (fl). (8.13) 
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Problem 2. Find a function v such that 

v £ V, (8.14) 
Y(h)-Y(v) + (L 2 (v),h-v)>(K + F,h-v), heV. (8.15) 

Here, the operator : X — ► X* is given by 



<9u 



i 



(u),h)=2 / e(|i?|,x)(A + + u)) 2 £ij{u + u)Eij(h)dx, u,heX. 



n 



(8.16) 



If (va, Pa) is a solution of Problem 1, then the pair (u + v\, p\) is a generalized solution of 
(4.1)-(4.4) with the viscosity function defined by (8.2). If v is a solution of Problem 2, then 
u + f is a generalized solution of (4.1)-(4.4) with the viscosity function defined by (8.1). 

Theorem 8.1. Suppose that the conditions (4.6), (4.7), (C4a) are satisfied, and let the 
function ipi satisfy both (CO) and one of the conditions (Cla), (C2a), (C3a). Then, for an 
arbitrary A > there exists a unique solution [v\, p\) of (8.11) -(8.13). Moreover, there 
exists a unique solution v of (8.14) -(8.15). In addition v\ — ^ v in V as A — > 0. 

Proof The existence of a solution (v\, p\) of (8.11)-(8.13) for an arbitrary A > follows 
from Theorem 5.1. Let (v{, p\) and {v\, p\) be two solutions of (8.11)-(8.13). By (8.12) 
we obtain 

{^(vl) + LM)-^(vl)-L 2 (vl)y x -vl) =0. (8.17) 



By Lemmas 3.1-3.4, the operator ^^-+L 2 is strictly monotone . Consequently (8.17) implies 
v\ = v\, whence p\ = p\. 

By Theorem 7.1 from the sequence {v\} a subsequence, again denoted by {v\}, can be 
selected such that v\ — 1 v in X where v is a solution of (8.14), (8.15). Also, Remark 7.1 
infers 

veV, Y(h) + $i(h)-Y(v)-$i(v)>(K + F,h-v), heV, (8.18) 

where 

^ U) = JnUo M^\E\,x)d^jdx, ueV, (8.19) 



and 



(^l(u),h)=(L 2 (u),h). (8.20) 



In addition, the problems (8.14), (8.15) and (8.18) are equivalent. The functional Y is 
convex, and the functional <E>i is strictly convex. Therefore the functional 

y 1 (u) = Y(u) + <!> 1 (u)-(K + F,u), ueV, 

is strictly convex, and if v\, v 2 are two solutions of the problem (8.18), then we have 

Z Z Z h£V 

Whence v\ = v 2 . ■ 
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9. Numerical solution of stationary problems. 

9.1. Theorem on convergence. Let {X m }, {N m } be sequences of finite-dimensional sub- 
spaces in X and L 2 (fi), respectively, such that (3.38), (3.39), (3.43) and (5.11) hold. We 
denote 

T dY x 



du 



+ L 2 , (9.1) 



where L 2 and are the operators defined by (8.10), (8.16). 

We seek an approximate solution (v m , p m ) of problem (8.11)-(8.13) of the form 

( v m,Pm) S X m X N m , (9.2) 

(L 3 (v m ),h)-(B mPm ,h) = (K + F,h), heX m , (9.3) 
(B m v m ,q)=0, qeN m . (9.4) 

We remind that the operator B m is defined by (3.40) and B m is the adjoint operator of B m . 

Theorem 9.1. Suppose that conditions (4.6), (4.7), (C4a) are satisfied, and let the function 
fa satisfy both (CO) and (C3a). Let also (3.38), (3.39), (3.43) and (5.11) are fulfilled. Then, 
for an arbitrary m £ N there exists a unique solution of (9.2) -(9.4) and 

v m v x in X, p m p x in L 2 (fi), (9.5) 

where v\, p\ is the solution of (8.11)-(8.13). If in addition ipi satisfies (Cla) or (C2a), 
then 

v m -> v x in X, (9.6) 

Pm ~^Px in L 2 (n). (9.7) 

Proof. The existence of a unique solution (v m ,p m ) of the problem (9.2)-(9.4) and the 
relations (9.5) follows from Theorems 5.1, 8.1. The following equalities also arise from the 
proof of Theorem 5.1 (see (5.35), (5.36)) 

(L 3 (v m ),v m ) = (K + F,v m ) -► (K + F,v x ), (9.8) 
lim(L 3 (v m ),h)-(B*px,h) = (K + F,h), heX. (9.9) 

(8.12), (9.8), (9.9) yield 

lim(L 3 (t> m ) - L(v x ),v m -v x ) = 0. (9.10) 

Assume that (C2a) is satisfied; if ipi meets (Cla) it meets also (C2a) (see Subsection 2.2). 
Then observing Lemmas 3.2, 3.4, we obtain 

{L 3 (v m ) - L 3 (v x ),v m - v x ) > 2a 3 / [I(u + v m )^ - I{u + v x )^} 2 dx, 

and (9.10) implies 

/ [I(u + v m )^ ~ I(u + vx)^] 2 dx 0. (9.11) 
Jn 

From here, taking into account that the function u —> JqU 2 dx is a continuous mapping 
from L 2 (f2 into R, we obtain 

/ I(u + v m )dx -> / I(u + v x )dx. (9.12) 
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It is obvious that 



\v m - V\\\x 



n IV n IV 

= / 52 ( £ ij( v m -v\)fdx = \ hj(u + v m ) - eij(u + v x )] 2 dx 

I Jn i,j=l 

= / I(u + v m )dx + I I(u + v\)dx — 2 / £jj(u + v m ) Eij{u + v\)dx. 

Jn Jn Jn (g 

By virtue of (9.5), (9.12) the right-hand side of (9.13) tends to zero. Therefore, (9.6) holds 
true. 

It follows from (8.12) and (9.1) that 

(L 3 (v x ),h)-(B*p x ,h) = (K + F,h), heX m , (9.14) 

By (9.3) and (9.14) we get 

(B*( Prn -n),h) = (L 3 (v m )-L 3 (vx),h) + (B*(p x -n),h), h £ X m , \i £ N m . 

This equality, together with (3.43), yields 

„ „ , (B*(p m - n),h) 

Pm ~ H l 2 (Q) < sup — — 

hex m 0\\h\\x 



where 



Hence 



<f3~ l {\\Lz{vm)-Lz{vx)\\x* +4px-^\\l 2 (Q)), ^4, (9.15) 
c = \\ B *\\c(L 2 (n),x*) = \\ B \\c(x,L 2 (n))- 

\\P\ -Pm\\L 2 (n) < \\P\ ~ ^\\l 2 (Q) + \\Pm ~ HU 2 (ft) 

< P^WL^Vm) - L 3 (v x )\\x* + (c/T 1 + l) inf \\p x - /i|| La(n ). (9.16) 

Lemmas 3.2, 3.4 and (9.6) imply L 3 (v m ) — > £3(1^) in X*, and (9.7) follows from (9.16) and 
(3.39). 

9.2. A saddle-point approach. We introduce two functionals 
J : X -> R, * : X x £ 2 (fi) -»• R defined by 

■/(«)= / ( / [e(| J E|,x)(A + 0"5 +Vi(e,I^U)]^)dx-(^ + F,n), 

Wo ' ( 9 . 17 ) 

1/) = J(u) - {B*v, u). (9.18) 

Problem: find a saddle-point of the Lagrangian ^, i.e. 

v x ,px£X xL 2 (Q), (9.19) 

¥(v A ,i/) < *(«a,pa) < *(«,pa), ^g^, (9.20) 

Theorem 9.2. Suppose that conditions (4.6), (4.7), (C4a) are satisfied, and let the function 
ipi meet (CO) and one of the conditions (CTa), (C2a), (C3a). TTten, for an arbitrary A > 
there exists a unique solution vx, Px of problem (9.19), (9.20) which is the solution of problem 
(8.11)-(8.13), and the problems (8.11)- (8.13) and (9.19), (9.20) are equivalent. 
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Proof. Let h be an arbitrary fixed element of X. Consider the function g : t — > g(t) = 
+ th,p\), t E K. It follows from the second inequality in (9.20) that g(0) < g(t), 
teR. The function g is differentiable in R, and so (0) = 0. This equality is equivalent to 

(8. 12) . By Lemmas 3.1-3.4 the operator L 3 is monotone, and so the functional u — > ^(n, g) 
is convex for an arbitrarily fixed g G L 2 {p) (see e.g. [5]). Therefore, the minimum of the 
functional u — > \l/(u,p,\) is characterized by (8.12). 

The first inequality in (9.20) gives (Bv\,p\ — v) < 0, for all v G L2(fi), and so we get 

(8.13) . However in this case ^(v\,p\) = $f(v\,v) for all v G Z,2(fi). Thus, the problems 
(8.11)-(8.13) and (9.19), (9.20) are equivalent. The existence and uniqueness of the solution 
of (9.19), (9.20) follows from Theorem 8.1. ■ 

Now we introduce the augmented Lagrangian as follows: 

*!(u,i/) = *(u,i/) + T -{Bu,Bu), (u,u) G X x L 2 (n). (9.21) 

It is obvious that the pair (v\,p\) is also the saddle point of the functional ^1, where r is 
an arbitrary positive constant. 
We study 

Algorithm of the augmented Lagrangian: find a sequence {v m ,p m } satisfying 



(v m+1 ,p rn+1 ) G X x L 2 (fi), (9.22) 

(L 3 (v m+1 ),h)-(B*p m ,h)+r(Bv m+1 ,Bh) = (K + F,h), heX, (9.23) 

(Pm+1-Pm,v) + Pm{Bv m+1 ,v) =0, U = L 2 (Q). (9.24) 

Here p m is a positive constant. 



Theorem 9.3. Suppose that conditions (4.6), (4.7), (C4a) are satisfied, and let the function 
ipi meets (CO) and one o/ i/ie conditions (Cla), (C2a). Assume that (vo,po) is an arbi- 
trary pair in X x £2(0). Then, /or an arbitrary m there exists a unique pair (v m+ i,p m+ i) 
satisfying (9.22) -(9.24). Moreover, if 

< inf p m < sup p m < 2r, (9.25) 

then 

v m -»■ u A in A", p m -> p A in L 2 (fi), (9-26) 
where (v\,p\) is the solution of (8.11) -(8.13). 
Proof. We set 

u m = v m - v x , q m =p m - Pa- (9-27) 

By subtracting (8.12) from (9.23) and (8.13) from (9.24), we get 

(L 3 (v m+1 )-L 3 (v x ),h)+r(Bu m+1 ,Bh) = (B*q m ,h), h G X, (9.28) 
{q m +i - q m ,v) = -p m (Bu m+1 ,u), v G L 2 (0,). (9.29) 

Taking v = 2g m+ i in (9.29), we obtain 

lkm+l||| 2 (Q) — lkm|li 2 (n) + lkm+1 — 9m||| 2 (Q) = ~ 2pm(Bu m+ i, q m +l)- 

(9.30) 
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Take h = u m+1 in (9.28). Then (9.28)-(9.30) give 



II II 2 _ II II 2 _i_ II _ II 2 

\\Qm+l\\L 2 (n) WQmW L 2 (fy + WQm+l <lm\\L 2 (Cl) 

+2p m (L 3 (v m+1 ) - (L 3 (vx),u m+1 ) + 2p m r\\Bu m+ i\\\ 2{u) 

= -2p m (Bu m+1 ,q m+1 - q m ). (9.31) 



(9.31) and Lemma 3.2 (see (3.15)) imply 

ll ll 2 ll ll 2 j_ ll ll 2 

\\1m+l\\L 2 (Q) ~ \\1m\\L 2 (Ci) + Wm+1 ~ Qm\\ L2 (fl) 

+^p m a 3 / [I(u + v m+1 )i - I(u + v\)^] 2 dx 
Jn 

+2p m r\\Bu m+1 \\l 2(n) < 2p m \\Bu m+1 \\ L2 ^)hm+i ~ q m \\L 2 (n)- (9-32) 
By applying the inequality 2ab < a 2 + b 2 to the right-hand side of (9.32), we obtain 

lkm+i||l 2 (Q) - lkm||| 2 (Q) +4p m a 3 J [I(u + v m+1 )? - I(u + v x )^] 2 dx 

+ Pm (2r - p m )\\Bu m+1 \\l 2{n) < 0. (9.33) 

By virtue of (9.25) there exists 5 > such that p m (2r — p m ) > 6 for any m, and by (9.33), 
\\q m \\L 2 (n) > \\Qm+i\\L 2 (n)- Thu s, the sequence {||<?m||| 2 ( Q) } converges, i.e. lim ||<?m||| 2( n) = 
a > 0, and (9.33) yields 

/ [I(u + t> m+1 )5 - I(u + vx) 1 *} 2 dx 0, (9.34) 
Bu m -► in L 2 (n). (9.35) 

Since the function u —> JqU 2 dx is a continuous mapping from L,2(£l) into M we obtain from 
(9.34) that 

/ I(u + v m )dx — > / I(u + v\)dx. (9.36) 
Jn Jn 

Therefore 

\\v m \\x < c, m e N, (9.37) 

and by (9.23), (3.37) we get ||pm||L 2 (n) < c f° r an m - Therefore, a subsequence {f^,^} can 
be extracted such that v v — ^ t>o in X, p T; — ^ po hi ^2(^)- We pass to the limit by analogy 
with the above. Then we get vo = v\, po = p\. Due to (9.36) and by the uniqueness of the 
solution of (8.11)-(8.13), we obtain by analogy with the above (see (9.13)) that 



v x in X. (9.38) 

It follows from (9.28) and (3.37) that 

\\L 3 (v m+ i) - L 3 (v x ) +rB*Bu m+1 \\ x * = \\B*q m \\ x * > Pi\\Qm\\L 2 (0)- 

This inequality, together with (9.35) and (9.38), yields q m — > in L,2{£1). Therefore (9.26) 
holds true. ■ 
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9.3. Solving a nonlinear problem. We consider two methods for solving the nonlinear 
problem (9.23), namely the Birger-Kachanov method and the contraction method. Both 
methods transform a nonlinear problem into a sequence of linear problems. 
We consider the problem: find a function u satisfying 

ueX, (L 3 (u),h) +r(Bu,Bh) = (f,h), V/i G X, (9.39) 

where / G X*. 

For an arbitrary v e X we define the operator M(v) G £{X, X*) as follows: 

(M(v)w,h) = 2 / [e(\E\,x)(X + I(u + v))~h ij (u + w)e ij (h) 
Jn 

+tpi(I(u + v),\E\,x) Sij(u + w) eij{h)]dx + r{Bw, Bh). (9.40) 

The Birger-Kachanov method consists in constructing a sequence {u m } such that 

u m+l G X, (M(u m )u m+l , h) = (/, h), V/i G X. (9.41) 

The conditions for the convergence of the Birger-Kachanov method in the general situation 
were established in [4]. 

From the known results [4], [13], [16], the next theorem follows. 

Theorem 9.4. Suppose the conditions (4.6), (4.7), (C4a) are satisfied. Assume that tpi 
is a nonincreasing function meeting the conditions (CO) and (Cla). Let also A > an uo 
be an arbitrary element of X. Then for any m there exists a unique solution v m+ \ of the 
problem (9.41) and u m — > u in X , where u is the solution of (9.39). 

Consider now the contraction method. Let A be a linear continuous selfadjoint and 
coercive mapping from X into X* , i.e. 

(Au,h) = (u, Ah), \(Au,h)\ <bi\\u\\x\\h\\x, u,h£X, 

(Au,u) >b 2 \\u\\ 2 x , (9.42) 
where b\, 62 are positive constants. By (9.42) the expression 

|H|i = (Au,u)^ (9.43) 

defines that norm in X that is equivalent to the norm ||.||x and to the norm of H 1 ^) 11 . By 
X\ we denote the space X equipped with the scalar product 

(«, h) Xl = (Au, h) (9.44) 

and with the norm (9.43). 

We study the following iterative method: 

u m+i G X, (Au m+1 ,h) = (Au m ,h) -t[L 3 (u m ,h) +r(Bu m ,Bh) - {f,h)], V/i G X, 

(9.45) 

where t is a positive constant. 

We may define the operator A by 

(Au,h) = J^g^^dx, u,heX, (9.46) 

where g G C(f2), g(x) > cq > for all x6!!. Then, taking h = (hi, 0) and h = (0, h 2 ) in the 
case that $7 C R 2 , and h = (hi, 0,0), h = (0,/i 2 ,0), h = (0,0, h 3 ) for Q C M 3 , where hi are 
arbitrary functions from [/ = {77 G if 1 (f2), 77 L = 0}, we split the problem (9.45) and obtain 
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independent problems for calculation u m i, i = 1, • • • , n. Such a split is very convenient for 
computations. 

Lemma 9.1. Suppose the conditions (4.6), (CO), (Cla), (C4a), (9.42) are satisfied. Then 

(L 3 (v) - L 3 (h),v- h) > gi \\v - h\\l (9.47) 
\\L 3 (v) - L 3 (h)\\ x * < q2\\v - h\\!, (9.48) 

where 

qi = Atifer 1 ' 

q 2 = {H2 +4a 5 A-i)6^. (9.49) 

Hi, fi2 are defined by (3.8). 

Proof. Taking into account the inequalities of (9.42), and applying Lemmas 3.1 and 3.4, 
we obtain (9.47), (9.48) with qi , q 2 defined by (9.49). 

Theorem 9.5. Suppose the conditions (4.6), (CO), (Cla), (C4a) are satisfied. Let f G X* 
and the operator A £ £(X, X*) meets (9.42). Let also A > and uq be an arbitrary element 
ofX. Then for t G (0, 2qiq^ 2 ), where 

q 3 = q 2 +r\\B*B\\ c{XltXl) , (9.50) 

and for any m there exists a unique solution u m+ \ of problem (9.45) and the following 
estimate holds 

k(t) m 

IK-Hli < l _ k{t) \\L 3 (u ) + rB*Buo-f\\x*, (9.51) 

where 

k(t) = (1 - 2 qi t + qlt 2 )^ < 1, (9.52) 

and u is the solution of (9.39). 

The function k takes its minimal value k(to) = (1 — q\ q 3 2 )z at the point to = q\ q% 2 . 

Proof. Let N = L 3 + rB*B. By Lemma 9.1 we have 

(N(v) - N(h), h-v)> qi \\v - h\\\, 

\\N{v)-N{h)\\xt<qz\\v-h\\ u v,heX, (9.53) 

where g 3 = q 2 + r\\B*B\\ ciXuX * ) . 

Denote by J the Riesz operator J 6 C(Xf,Xi) that is defined as follows 

(Jg,h) Xl = (g,h), geXl h G X\, \\ Jg^ = \\g\\ x . . (9.54) 

It is obvious that the problem (9.39) is equivalent to finding a fixed point u = Ut(u), where 
Ut-.X^Xi, 

U t (h) =h-tJ(N(h) - /). (9.55) 

By (9.53)-(9.55) we have 

\\U t (v) - U t (h)\\ 2 = \\v-h- U(N(v) - N(h))\\ 2 
= \\v - h\\l - 2t(N{v) - N(h), v-h) + t 2 \\N(v) - N(h)\\ 2 x * 
< \\v-h\\ 2 -2t qi \\v-h\\ 2 +t 2 q 2 \\v-h\\ 2 = (k^Wv-hh) 2 , 
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where k(t) is defined by (9.52) and k(t) < 1, if t € (0, 2gi(/^ 2 ). Therefore the mapping Ut 
is a contraction, and the existence of a unique solution u of problem (9.39) and estimate 
(9.51) follow from the fixed point theorem (see e.q. [21]). ■. 

Remark 9.1. Equalities (9.49) and (9.50) imply that q^ — > oo as A — > 0. Therefore 
k(to) = (1 — qfq^ 2 )^ the minimal value of k(t) tends to unit as A — > 0, and the iterative 
method (9.45) provides slow convergence at small value of A. The reason of this is that 
the differentiable functional Y\ tends to nondifferentiable functional Y (see (8.8), (8.9)) as 
A -► 0. 

9.4. Solving the problem (5.8)— (5.10). For the case that the conditions of Theorem 5.2 
are satisfied, the operators and L are strictly monotone and Lipschitz continuous. There- 
fore the algorithm of the augmented Lagrangian (see (9.22)-(9.24)) can be used for the 
solution of the problem (5.8)-(5.10), and the corresponding nonlinear systems can be solved 
by the Birger-Kachanov method and by the contraction method. 

Let us consider the general case that the conditions of Theorem 5.2 are not satisfied. 
Let {Xi}*=i m \ {^i}f=i"^ be bases in the spaces X m and N m , respectively. Let also k = 
ki(m) + fc 2 (m). Define a mapping M : R h -> R h as follows: 

R k 3c = { Cl } k t=l -^M(c) = {Mi(c)}U, (9-56) 

where 

fci(m) fei(m) fei(m) 

j=i j=i j=i 

k 

"(C c 3Vj-ki{m),Xi)-(K + F,Xi), i = l,...,fci(m). (9.57) 

j=fci(m)+l 

fei(m) 

Mi(c) = (#m Y c j'Xj,%-fci(m)), i = fci(m) + 1, (9.58) 
j'=i 

It is obvious that the problem (5.8)-(5.10) is equivalent to the following one: find c = 
(6i, . . . , Cfc) such that 

c£K fc , M(c) = 0. (9.59) 

Define the functional 

k 

$ 2 (c) = ^(M(c)) 2 , c G (9.60) 
i=i 

The problem (9.59) is equivalent to the following one: 

c G M fe , $ 2 (c) = min $ 2 (c) = 0. (9.61) 

In the case that the functions tp an d b are continuously differentiable the functional $ 2 is 
continuously differentiable in R fc , and gradient method can be applied for calculation of a 
solution of the problem (9.61). Derivative of the mapping v — > ^^(v,v) + L(i>) is defined 
in 5.2 (see (5.43)-(5.48)). 
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10. Stationary problem with consideration for the inertia forces. 

nonhomogeneous problem. 

10.1. Basic equations and auxiliary results. The equations of motion with regard for 
the inertia forces read as follows: 

puj + ^ 2 —[(p(I(u), \E\, fi(u, E))sij (u)] = Ki in tt, i = l,...,n. 

OXj OXi OXj (10.1) 

The condition of incompressibility is 

divu = 0. (10.2) 
We assume that velocities are specified on the boundary S of Q, i.e. 

u = u. (10.3) 

s 

We assume also 

(C7): $7 is a bounded domain in R n , n = 2 or 3. The boundary S of f2 belongs to the 
class C 2 and consists of / connected components F±, . . . , Fi (I > 1), 
and suppose that 

u£H^(S), UiVids = 0, i=l,..., I, (10.4) 
Jr, 

and the function ip is defined by (2.12). 

The following lemma follows from the known results (see e.g. [11, 23]). 

Lemma 10.1. Suppose that the conditions (CI) and (10.4) are satisfied. Then there exists 
a function u such that 

u = curl r/, rj = (771, . . . ,f] n ) £ H 2 (9) n , u\ s = u, (10.5) 

moreover for an arbitrary a > one can choose a vector-valued function n such that 

We set 



UiVj\\ L2 (si)<a\\v\\ x , veHn(n) n , i,j = l,...,n. (10.6) 



q(u,v,w) = p / Uk-^-Widx, u,v,w £ H 1 (fl) n . (10-7) 
Jn oxk 

Obviously 

n 

\q(u,v,w)\ < p \\ u kh 4 (Q) g^- T ^ m \\wi\\ Li (n), (10.8) 



i,k=l 

Therefore, the trilinear form q is continuous in i7 1 (0) n x H 1 (Q) n x H 1 ^) 71 . 
We consider the following spaces 

X = H^(n) n with the norm || ■ \\ x = \\ ■ \\ x , (10.9) 

V = {w £ X, &xvw = 0} with the norm || • ||y = || • \\x, (10.10) 



M = {w £ L2(fi), / wdx = 0} with the norm || • \\n = \\ ■ ||L2(fi). 
Jn 

It is easy to verify that 



(10.11) 

q(z,w,h) = —q(z,h,w), z£V, w,h£H 1 (Q) n , n = 2 or 3, q(z,h,h)=0. 

(10.12) 
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Define a trilinear form q± as follows: 

qi(v, w, h) = — q(v, w, h) — — q(v , h, w), v, h, w G H 1 ^) 11 . (10.13) 
It is evident that 

qi(v,h,h) =0, v^eH 1 ^)", (10.14) 

qi(v,w,h) = -qi(v,h,w), v,h,w £ H l ($l) n . (10.15) 
Lemma 10.2. Suppose that 

v k — 1 f in X, (10.16) 

v k -> v in L 4 (n) n . (10.17) 
Then, for an arbitrary fixed h G iJg(n) n i/ie following relations hold 

Qi(vk,Vk,h) -> qi(v,v,h), (10.18) 

qi(u,v k ,h) -> qi(u,v,h), (10.19) 

qi(v k ,u,h) ->■ qi(v,u,h). (10.20) 

in i/iis case ifv&V, then 

qi(v,v,h) = q(v,v,h), qi(u,v,h) = q(u,v,h), qi(v, u, h) = q(v, u, h). 

(10.21) 

Proof. We have 

- wjM 2 < ( _£(«fcj - ^o 4 dx Y{J n h i dx )^ ( io - 22 ) 

where v^- are the components of the vector-valued function v k . By (10.17) the left-hand 
side of (10.22) tends to zero. 
Therefore, 

Vkjhi — > fj/ii in L2(fi) as fc — > oo, i, j = 1, . . . ,n. (10.23) 
(10.16) and (10.23) yield 

q(v k ,v k ,h) -> q(v,v,h). (10.24) 

Application Green's formula gives 

q(vk,h,v k ) = p [ v kj ^-v ki dx = -p j hi-^-(v kj v ki )dx, 

and by analogy with the stated above we obtain 

q(v k ,h,v k ) -> q(v,h,v). 

Therefore (10.18) holds. It is evident that (10.19), (10.20) follows from (10.16) and (10.17). 
In the special case that v G V we have 

/* /* dh' 

q(v,v,h) = p / Vj—-^-hidx = —p / Vj—-^Vidx = —q(v,h,v). (10.25) 
°Xj is] 

(10.13) and (10.25) imply (10.21). 

Lemma 10.3. Suppose the conditions (C7) and (10.4) are satisfied. Then for an arbitrary 
£ > one can choose a vector valued function n swc/t i/tai (10.5) is satisfied and in addition 

\qi(v,u,v)\<av\\ 2 x, veX. (10.26) 
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Proof. By Green's formula we obtain 
qi(v,u,v) = \p j 



dui 



dvi 



n 

UiVi 



3 dx 3 



u. 



dx 



dx. 



v £ X. 



It follows from here that 



\qi(v,u,v)\ <ci||w||x 11^11^2(0); 

and Lemma 10.3 follows from this inequality and Lemmas 10.1. 

10.2. Boundary value problem. We consider the problem: find v satisfying 
v e V 



(10.27) 



\^r( v > v )> w ) + ( L ( v )> w ) + v, w) + q(u, v, w) + q(v, u, w) - (y, w) = 0, w £ V, 
V ° n ' (10.28) 



where 



(y, w) = KiWi dx - 

Jn Jn 



Uj— — Wi dx. 



(10.29) 



It follows from (4.7) and (10.5) that y G X* . The left-hand side of (10.28) belongs to the 
polar set 

V° = (f,w) =0, WGV}. 

Therefore, there exists a function p G N such that the pair (v, p) is a solution of the following 
problem: 

(v,p)eXxM, (10.30) 

'dJ x 



+q(v, u, w) - (B*p, w) = (y,w), we X, 
(Bv,j)=0, 7GAA. 



(10.31) 
(10.32) 



By use of Green's formula it can be seen that, if (v,p) is a solution of problem (10.30)- 
(10.32), then (u,p) with u = u + v is a solution of problem (10.1)-(10.3) in the sense of 
distributions. On the contrary, if (u,p) is a classical solution of problem (10.1)-(10.3), then 
the pair (v,p) with v = u — u is a solution of problem (10.30)-(10.32). 

Let {X m }, {M m } be sequences of finite-dimensional subspaces in X and N which satisfy 
the following conditions 



lim inf \\w — z\\x = 0, 

m— «x> zeX m 

lim inf \\h - y\\ L2 (n) = 0, 

m-»oo yeNm 



w G X, 
h e M, 



inf sup 



(B m w,n) 



veNm wex m \\ w \\xw\\L 2 (n) 

X m C X m +1, M m C Alm+l 



> (3 > 0, men, 



(10.33) 
(10.34) 

(10.35) 
(10.36) 
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We introduce the spaces V m and \>° m by 

V m = {ueX m , (B m u, 7 )=0, 7 eAA m }, (10.37) 

V; = {?G^, (9,u)=0, nGV m }. (10.38) 
Define approximate solutions of problem (10.30)-(10.32) as follows: 

(v m ,p m ) £ X m xM m , (10.39) 

d J 

(-^("m.fm),^ + (L(u m ),w) + <7i(v m ,v m ,w) +gi(tt,w m ,u;) 

+gi(f m ,n,w) - (B^p m ,w) = (y,w), w e X m , (10.40) 
{B m , v m , 7 ) = 0, 7 € M m . (10.41) 

Theorem 10.1. Suppose that K 6 L2(^) n and i/ie condition (C4) zs satisfied. Let the 
function^ meets one of conditions (CI), (C2), (C3) (V replaced byip). Assume that (C7), 
(10.4), (10.5) and (10.26) with £ < ai are fulfilled. Let also {X m } and {J\f m } be sequences 
of finite- dimensional subspaces in X and N respectively, such that (10.33) -(10.36) hold. 
Then, for an arbitrary fixed A > and an arbitrary m 6 N there exists a solution (v m ,p m ) 
of problem (10.39)-(10.41) ; and a subsequence {vk,Pk} can be extracted from the sequence 
{v m ,Pm} such that Vj~ — v in X, pk — 1 p in N , where v,p is a solution of the problem 
(10.30) -(10.32). 

Proof. Define an operator M : X — > X* by 
(M(g),w) = (-^-(g,g),w) + ( L (g),w) + qi{g,q,w) +q 1 (u,g,w) +qi(g,u,w), g,w £ X. 

\ an j (\§A1) 

It follows from (10.39)-(10.42) that v m is a solution of the problem 

v m eV m , {M(v m ),w) = (y,w), weV m . (10.43) 

Bearing in mind (10.14) and that £ < a± in (10.26) we obtain by analogy with the proof of 
Theorem 5.1 (see (5.13), (5.15), that 

z(e) = (M (e), e) - {y, e) > ai||e|& - c||e||x, e€X, A > 0. (10.44) 

Therefore, z(e) > for ||e||x >r=-^, and there exists a solution of (10.43) with 

\\v m \\x<r, \\M(v m )\\x* <ci m G N. (10.45) 

For an arbitrary f <E X* we denote by Gf the restriction of / to X m . Then Gf G X^, and 
by (10.43) we obtain 

G(M(v m ) -y)e V° m . (10.46) 

Therefore, there exists a unique p m € N m such that (10.40) is satisfied (see Lemma 3.6), 
and (10.35) yields 

\\Pm\W < C2- (10.47) 
By (10.45), (10.47) we can extract a subsequence {vk,Pk} such that 

u fc v in AT, (10.48) 

Vk — ► f in L 4 (f2) n and a.e. in 17, (10.49) 

M(vk) ^ f] in X*, (10.50) 

p fc p in AT. (10.51) 
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Next we use Lemma 10.2 and by analogy with the proof of Theorem 5.1 we pass to the limit 
in (10.40), (10.41) with m replaced by k, and obtain 

(M(v),w)-(B*p,w) = (y,w), weX. (10.52) 
(Bv,j) = 0, 7 gAA. (10.53) 

Taking into consideration that (10.21) is satisfied, for v G V, we get that pair (v,p) is a 
solution of the problem (10.30)-(10.32). 

11. Stationary problem with consideration for the inertia forces. Mixed 

PROBLEM. 

11.1. Formulation of the problem and an existence result. As before we consider that 
Si and S 2 are open subsets of the boundary S of f2 such that Si is non-empty, Si n S2 = 
and Si U S 2 = S. We study the problem on searching for a pair of functions (u,p) which 
satisfy the motion equations (10.1), the condition of incompressibility (10.2) and the mixed 
boundary conditions, wherein velocities are specified on Si and surface forces are given on 
S 2 , i.e. 



u 



= u, (11.1) 

01 



[-p5 ij + 2tp{I{u),\E\,n(u,E))e ij {u)}v j = F u i = l,...,n. (11.2) 

S2 

It is obvious that in the special case that S 2 is an empty set this problem transforms into 
the problem considered in Section 10. 

We assume that (p is defined by (2.12) and u G Hz {Si). Then there exists a function u 
satisfying (4.6). 

Let us define operators Mi : X — > X*, M 2 : X — > X* and an element x G A* as follows: 

(MiH,s) = (^( w ,w),gj +{L{w),g), w,g G X, (11.3) 
{M 2 {w),g) = q{u,w,g) +q{w,u,g), w,g£X; (11.4) 

(X,ff)= / Kigidx+ Figids - q{u,u,g), g £ X. (11.5) 

We consider the problem 

M^x^), (11-6) 

(Mi(u),«;) + (M2(i;),«;) + g(i;,v,«;)-(S*p,w) = (x.«'), w G X, (11.7) 

(Bu,7) = 0, (11.8) 

By using Green's formula, one may show that, if {v,p) is a solution of the problem (11.6)- 
(11.8), then {u,p) with u = u + v is a solution of the problem (10.1), (10.2), (11.1), (11.2) 
in the distribution sense. On the contrary, if (u,p) is a solution of (10.1), (10.2), (11.1), 
(11.2) such that (11.6) holds with v = u — u, then {v,p) is a solution of the problem (11.6) 
-(11.8). Define also the following constants 

ri = sup q{w, w, w), 
w£V,\\w\\x<l 

r 2 = inf (M 2 {w),w), 
w£V,\\w\\x<l 

r 3 = sup \(x,w)\- (11-9) 

uiGV,||ui||x<1 
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Consider the space 

P = {w£H 1 (Q) n , dww = 0}. (11.10) 

The space P is presented in the form P = V © F -1- , where V is given by (3.2) and is the 
orthogonal complement of V in P. Evidently that the constants r 2 and r% in (11.9) depend 
on the function u satisfying (4.6). Let u\ be a function from V 1 - that satisfies (4.6). Then 
the function u = u\ + U2, where U2 is an arbitrary element of V, meets (4.6). We assume 
that there exists a function U2 £ V such that the following inequalities hold 

r 4 = ai + y > 0, rf > nr 3 , (11.11) 

where ai is the positive constants from (2.13). It is evident that (11.11) holds, if the norms 
of the functions K, F, and u are not large. 

Lemma 11.1. Suppose the condition (C4) is satisfied and the function ip meets one of the 
conditions (CI), (C2), (C3) (V replaced byip). Let also (11.11) holds. 
Then the following inequality is valid: 

0(w) = (Mi (w),w) + (M 2 (w),w) +q(w,w,w) - (x,w) > 0, 



(11.13) 



•/ r- \r J II II r 4 - y/r'i-rir 3 m ,„> 

if w £ V and \\w \\x = 1 — = £■ (11.12) 

n 

Proof. It follows from (11.9) and (11.12) that 

/?(w) > /?i(||w||x) = (2ai + r 2 )|H|x - - r 3 ||u;||x, w £ V. 

Consider the quadratic equation 

2r 4 y - ny 2 - r 3 = 0, 

r 4 being defined in (11.11). Its roots are those of the equation 0i(y) = and they are equal 
to 

?~4 - V r 4~ r l r 3 r 4 + \/ r 4 ~ r l r 3 

yi = , V2 = • 

n n 

If (11.11) holds, then y\ and y2 are real and (3\(y) > for y £ [yi,y2\- Therefore (11.11) 
yields (11.12), and the lemma is proved. ■ 

Theorem 11.1. Let Q be a bounded domain in M. n , n = 2 or 3 with a Lipschitz continuous 
boundary S. Suppose the condition (C4) is satisfied and the function ip meets one of the 
conditions (CI), (C2) ; (C3) (ip replaced byip). Let also (11.11) holds. Then, for an arbitrary 
A > there exists a solution of the problem (11.6) -(11.8). 

Proof. It follows from (11.6)-(11.8) that the function v is a solution of the problem 

v £ V, 

(M 1 (v),w) + (M 2 (v),w)+q(v,v,w)=(x,w), w £ V. (H-14) 
Let {V m } be a sequence of finite-dimensional subspaces of V such that 

lim inf ||«; — z\\x = 0, w £ V, (11.15) 

m^oo zeV m 

V m C V m+1 . (11.16) 
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We search for the Galerkin approximations v k satisfying 

(Mi(v m ),w) + (M 2 (v m ),w) + q(v m ,v m ,w) = {x,w), w G V m . 

(11.17) 

By virtue of Lemma 11.1 there exists a solution of the problem (11.17) and ||w m ||x < £■ 
Thus we can extract a subsequence {v k } such that — 1 vq in V. We pass to the limit as 
k — > oo in (11.17), with m changed by fc. In this case by analogy with the stated above, (see 
the proofs of Theorems 5.1, 10.1 and Lemma 10.2), we use the methods of monotonicity and 
compactness. Thus, we get that the function v = vq is a solution of the problem (11.14). 
Now from Lemma 3.5 it follows that there exists a function p G L2(fi) such that the pair 
(v,p) is a solution of the problem (11.6)— (11.8). 

11.2. Approximation of the problem (11.6)— (11.8). Let {X m }, {N m } be sequences of 
finite dimensional subspaces of X and L2(fi) which satisfy the conditions (3.38), (3.39), 
(3.43) and (5.11). We search for an approximate solutions of the problem (11.6)-(11.8) in 
the form 

(v m ,p m ) G X m x N m , (11.18) 
(Mi(v m ),w) + (M 2 {v m ),w) + q(v m ,v m ,w) - (B^p m ,w) = (x,w), w G X m , 

(11.19) 

7 ) = 0, 7 G AT m . (11.20) 

From the point of view of applications, in particular, of computation, the problem (11.18)- 
(11.20) is considerably more preferable than (11.17). So, we study the question of conver- 
gence of the approximations {v m ,p m }. 
Define constants rji — 774 by 

771 = sup max \q(w,w,w)\, 
7/2 = inf min (M 2 (w),w), 

m w£d m 

7/3 = sup max \(x,w)\, 

m wGdm 

7/4 = 0! + ^, (11.21) 

where 

d m = {weX m , (B m w,j)=0, T^m, Hlx<l}. (H-22) 

Note that 7/1 > 77, t/ 2 < r 2 , 7/3 > r 3 (see (3.38), (3.39) and (11.9)). 
Define a mapping M 3 : X — > X* as follows: 

(M3(?j), 7D ) = (M 2 (t;), to) + g(u, u, ttj) = </(u, u, to) + q(v, u, w) + q(v, v,w), v, w G X. 

(11.23) 

Lemma 11.2. Suppose that (4.6) is satisfied and let 

{v k } C X, v k ^v in X. (11.24) 

Then 

lim(M 3 (v k ),w) = {M 3 (v),w), w £ X, (11.25) 
]im(M 3 (v k ),v k ) = (M 3 (v),v). (11.26) 
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Proof. It follows from (11.24) that 



Vk 



v in L^(Q), 



(11.27) 



and (11.25) arises from the proof of Lemma 10.2, see (10.24). 
We have 



+ 



q(vk,Vk,v k ) - q(v,v,v)\ 
dvki 



, dv ki 
(v k j - Vj)-^-Vkidx 



dv ■ \ f Qv ' 

-r- L )vkidx+ l Vj—^-^ki-v^dx. (11.28) 
in ^oxj oxj/ Jq oxj 

(11.27) implies (see (10.22)) that 

VjVki — ► VjVi in L 2 (fi). (11.29) 

By (11.24), (11.27), (11.29) each addend in the right-hand side of the equality (11.28) tends 
to zero. Therefore, 

limq(v k ,v k ,Vk) = q(v,v,v) 

and (11.26) holds true. ■ 

Theorem 11.2. Let Q be a bounded domain in W l , n = 2 or 3 with a Lipschitz continuous 
boundary S. Suppose the condition (C4) is satisfied and the function ip meets one of the 
conditions (CI), (C2), (C3) (<p replaced byip). Let also {l m },{JV m } be sequences of finite 
dimensional subspaces of X and £2(0) which satisfy the conditions (3.38), (3.39), (3.43) 
and (5.11). Finally, assume that 



774 > 0, r? 4 > 771773. 



(11.30) 



Then, for an arbitrary fixed A > 0, and for each m G N, there exists a solution of the 
problem (11.18)-(11.20), and a subsequence {(vk,Pk)} can be extracted from the sequence 
{(v m ,p m )}) such that Vk — v in X , pk — 1 p in L2{£1), where (v,p) is a solution of the 
problem (11.6) -(11.8). 

The proof of this theorem is analogous to that of Theorem 10.1, and it is not given because 
of this. 

Remark 11.1. In the general case we can consider that the density of an electrorheo- 
logical fluid depends on the module of the vector of electric field strength, i.e. p = p(\E\), 
and 



P2>p(y)>pi, ye 



(11.31) 



where pi, P2 are positive constants. It is easy to see that all results of Sections 10 and 11 
still stand valid in the case that p is a function satisfying the condition (11.31). 
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